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This  report  describes  the  operation,  maintenance  and  research  activities  at 
the  Norwegian  Seismic  Array  (NORSAR)  for  the  period  1  October  1982  to  31 
March  1983. 

The  uptime  of  the  NORSAR  online  detection  processor  system  has  averaged  95.7 
as  compared  to  92.3  for  the  previous  period.  Most  of  the  downtime  was  caused 
by  CPU  and  disk  drive  problems.  A  total  of  1208  events  were  reported  in 
this  period,  giving  a  daily  average  of  6.6  events.  The  number  of  reported 
events  per  month  varies  from  144  in  November  to  268  in  March.  There  have 
been  some  difficulties  with  the  communications  lines;  06C  was  down  at  the 
beginning  of  the  period  due  to  a  seriously  damaged  cable  which  was  finally 
repaired  3  November.  All  communications  circuits  were  turned  off  12  December 
in  connection  with  dismantling/removal  of  the  SPS.  There  were  also  some 
difficulties  with  the  telemetry  stations  at  02B  due  to  snow  and  ice  reducing 
the  charging  effect  of  the  solar  cells. 

The  DP  system  implemented  on  IBM  4331/MODCOMP  is  described  in  Section  III. 
Programs  and  operational  routines  used  on  the  SPS-system  have  been  converted 
Graphic  displays  ate  now  used  in  array  monitoring  routines.  A  data-loss 
problem  caused  by  sharing  of  disk  between  the  IBM  4331  and  4341  has  been 
corrected,  and  during  the  second  half  of  the  reporting  period,  the  shared- 
disk  principle  was  instigated  for  programs  as  well.  A  synchronization  problei 
in  the  MODCOMP  system  has  been  corrected  by  restarting  and  resynchronizing 
the  system  twice  a  week. 

Section  IV  describes  field  instrumentation  and  maintenance  activities  at 
the  NORSAR  Maintenance  Center. 

The  research  activity  is  briefly  described  in  Section  VI.  Subsection  1 
discusses  the  absorption  band  effect  on  long-period  body  waves.  Sub¬ 
section  2  discusses  the  present  seismic  evidence  for  a  boundary  layer 
at  the  base  of  the  mantle.  Subsection  3  presents  the  effect  of  aliasing  on 
the  NORSAR  detector  performance.  Subsection  4  presents  results  of  a  study 
of  seismic  noise  at  high  frequencies.  Subsection  5  discusses  power  spectral 
bias  sources  and  quantization  levels.  In  subsection  6  the  results  of  a 
study  of  magnitudes  from  P  coda  and  Lg  using  NORSAR  data  are  presented. 
Subsection  7  discusses  the  design  work  for  a  new  regional  array  in  Norway. 
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I.  SUMMARY 

This  report  describes  the  operation,  maintenance  and  research  activities  at 
the  Norwegian  Seismic  Array  (NORSAR)  for  the  period  1  October  1982  to  31 
March  1983. 

The  uptime  of  the  NORSAR  online  detection  processor  system  has  averaged  95.7%, 
as  compared  to  92.3  for  the  previous  period.  Most  of  the  downtime  was  caused 
by  CPU  and  disk  drive  problems.  A  total  of  1208  events  were  reported  in 
this  period,  giving  a  daily  average  of  6.6  events.  The  number  of  reported 
events  per  month  varies  from  144  in  November  to  268  in  March.  There  have 
been  some  difficulties  with  the  communications  lines;  06C  was  down  at  the 
beginning  of  the  period  due  to  a  seriously  damaged  cable  which  was  finally 
repaired  3  November.  All  communications  circuits  were  turned  off  12  December 
in  connection  with  dismantling/removal  of  the  SPS.  There  were  also  some 
difficulties  with  the  telemetry  stations  at  02B  due  to  snow  and  ice  reducing 
the  charging  effect  of  the  solar  cells. 

The  DP  system  implemented  on  IBM  4331/MODCOMP  is  described  in  Section  HI. 
Programs  and  operational  routines  used  on  the  SPS-system  have  been  converted. 
Graphic  displays  are  now  used  in  array  monitoring  routines.  A  data-loss 
problem  caused  by  sharing  of  disk  between  the  IBM  4331  and  4341  has  been 
corrected,  and  during  the  second  half  of  the  reporting  period,  the  shared- 
disk  principle  was  instigated  for  programs  as  well.  A  synchronization  problem 
in  the  MODCOMP  system  has  been  corrected  by  restarting  and  resynchronizing 
the  system  twice  a  week. 

Section  IV  describes  field  instrumentation  and  maintenance  activities  at 
the  NORSAR  Maintenance  Center. 

The  research  activity  is  briefly  described  in  Section  VI.  Subsection  1 
discusses  the  absorption  band  effect  on  long-period  body  waves.  Sub¬ 
section  2  discusses  the  present  seismic  evidence  for  a  boundary  layer 
at  the  base  of  the  mantle.  Subsection  3  presents  the  effect  of  aliasing  on 
the  NORSAR  detector  performance.  Subsection  4  presents  results  of  a  study 
of  seismic  noise  at  high  frequencies.  Subsection  5  discusses  power  spectral 
bias  sources  and  quantization  levels.  In  subsection  6  the  results  of  a 
study  of  magnitudes  from  P  coda  and  Lg  using  NORSAR  data  are  presented. 
Subsection  7  discusses  the  design  work  for  a  new  regional  array  in  Norway. 
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II.  OPERATION  OF  ALL  SYSTEMS 

II. 1  Detection  Processor  (DP)  Operation 

There  have  been  61  breaks  In  Che  otherwise  continuous  operation  of  the 
NORSAR  online  system  within  the  current  6-month  reporting  interval. 

Almost  half  of  the  downtime  in  this  period  was  caused  by  a  cut  of  the  main 
power  due  to  installation  of  a  new  transformer.  Most  of  the  stops  due  to 
DP  software  occurred  in  the  first  half  of  the  period;  the  last  month  had 
only  3  stops.  The  uptime  percentage  for  the  period  is  95.7  per  cent  as 
compared  to  92.3  for  the  previous  period. 

Fig.  II. 1.1  and  the  accompanying  Table  II. 1.1  both  show  the  daily  DP  downtime 
for  the  days  between  1  October  1982  and  31  March  1983.  The  monthly  recording 
times  and  percentages  are  given  in  Table  II. 1.2. 


The  breaks  can  be  grouped  as  follows : 

a)  Stops  related  to  possible  program  errors  51 

b)  Maintenance  stops  2 

c)  Power  breaks  6 

d)  Hardware  problems  2 


The  total  downtime  for  the  period  was  185  hours  and  50  minutes.  The  mean- 
time-between-failures  (MTBF)  was  2.9  days  as  compared  with  1.4  days  for  the 
previous  period. 


J.  Torstveit 


LIST  Cf  BREAKS  IN  CP  PRCCESSING  THE  LAST  HALF-YEAR 

DAY  START  STCP  CCPKENTS . . .  CAY  START  STOP  CCPHENTS 
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Table  II-l.l  List  of  breaks  In  DP  processing  In  the  period  1  Oct  82  -  31  Mar  83 


Month 

DP 

Uptime 

(hrs) 

DP 

Uptime 

<%) 

No.  of 

DP  Breaks 

No.  of 

Days  with 
Breaks 

DP 

MTBF* 

(days) 

Oct 

727.47 

97.8 

15 

12 

1.9 

Nov 

635.70 

87.7 

14 

13 

1.8 

Dec 

717.72 

95.8 

13 

9 

2.3 

Jan 

717.57 

96.4 

11 

8 

2.5 

Feb 

703.68 

97.6 

5 

4 

4.6 

Mar 

737.08 

99.1 

3 

3 

7.7 

4239.22 

95.7 

61 

49 

2.9 

*Mean-time-between-failures  *  (Total  uptime/No.  of  up  intervals) 

TABLE  II. 1.2 

Online  System  Performance 
1  October  1982  -  31  March  1983 


II. 2  Event  Processor  Operation 

In  Table  II. 2.1  some  monthly  statistics  of  the  Event  Processor  operation  are 
given: 


Teleseismic 

Core  Phases 

Sum 

Daily 

Oct  82 

152 

43 

195 

6.2 

Nov  82 

107 

37 

144 

4.8 

Dec  82 

124 

40 

164 

5.3 

Jan  83 

185 

33 

218 

7.0 

Feb  83 

183 

36 

219 

7.8 

Mar  83 

238 

30 

268 

8.6 

989 

219 

1208 

6.6 

TABLE  II. 2.1 


B.  Kr.  Hokland 
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II. 3  Array  Communication 

Table  II. 3.1  reflects  the  performance  of  the  communications  system  through¬ 
out  the  reporting  period. 

As  also  previously  stated,  the  table  not  only  reflects  irregularities  in 
the  communication  systems  itself,  but  in  addition  also  reflects  other 
prominent  conditions  such  as: 

CTV  power  failure 

Line  switching  during  program  debugging/corrections  (at  the  NDPC) 

-  Maintenance  visits 

NDPC-initiated  tests 
MODCOMP  resynchronization. 

Since  21  September  1982  the  MODCOMP  Communications  Processor  has  replaced 
the  SPS. 


Summary 

October :  06C,  which  had  been  down  since  week  38/82,  was  still  out  of 

operation  by  end  of  October  due  to  seriously  damaged  cable 

November :  06C  resumed  operation  (3  Nov).  A  few  resynch,  problems  observed 
in  connection  with  01B  in  the  beginning  of  November;  otherwise, 
reliable  system  operation. 

December :  12  Dec  all  communications  circuits  were  turned  off  in  connec¬ 
tion  with  SPS  dismantling/removal.  Due  to  power  line  damage 
in  the  02B  area  (11  Dec),  this  subarray  was  out  of  operation 
the  remaining  days  of  December  in  spite  of  strong  efforts 
from  the  local  power  plant's  crew. 

January :  Power  back  in  the  02B  area  (4  Jan),  but  due  to  SLEM  trouble, 
the  subarray  resumed  operation  5  Jan.  Found  spikes  in  the 
02C  data  (25  Jan),  and  this  initiated  a  test  and  investigation 
period  where  telestations  at  Llllestrem,  Hamar  and  Lillehamer 
were  involved. 


February  :  After  9  Feb  the  02C  system  was  back  In  operation. 

March:  Between  5  and  11  Mar  the  03C  communications  system  was  out 

of  operation  due  to  NTA  cable  work  in  the  Rena  area.  06C 
lost  synchronization  a  number  of  times  weeks  10-13. 

Table  II. 3. 2  indicates  distribution  of  outages  with  respect  to  the  individual 
subarrays. 

Miscellaneous 

Problems  occurred  in  connection  with  02B  telemetry  stations  11  Dec  due  to 
snow  problems  reducing  the  charging  effect  of  the  solar  cells  to  the 
batteries.  Operation  resumed  27  Dec.  The  same  problem  occurred  in  January 
(10-14),  but  this  time  only  telemetry  station  04  was  involved. 


► 


O.A.  Hansen 


OCT  (4) NOV  (4) DEC  (5) JAN  (4) FEB  (4)  MAR  (5)  AVERAGE  i  YEAR 

(4-31.10)  (1-28.11)  (29.11-2.1.83)  (3.1-30.1)  (31.1-27.2)  (28.2-3.4) 

>20  >200  >20  >200  >20  >200  >20  >200  >20  >200  S5o  7?7io  ^20  ^2nn 
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Week/ 

Subarray/per  cent 

outage 

Year 

01 A 

01B 

02B 

02C 

03C 

04  C 

06C 

40/83 

_ 

- 

— 

— 

_ 

100.0 

41 

- 

- 

- 

- 

0.3 

- 

100.0 

42 

- 

- 

- 

0.6 

- 

- 

100.0 

43 

- 

- 

- 

- 

- 

- 

100.0 

44 

- 

- 

- 

- 

- 

0.9 

1.2 

45 

- 

- 

- 

- 

- 

- 

0.6 

46 

7.1 

7.1 

7.1 

7.1 

7.1 

7.1 

6.5 

47 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

48 

49 

50 

- 

14.3 

- 

- 

- 

6.1 

- 

- 

100.0 

- 

- 

- 

- 

51 

- 

- 

100.0 

- 

- 

- 

- 

52 

- 

- 

100.0 

- 

- 

- 

- 

1/83 

- 

- 

39.3 

- 

- 

- 

- 

2 

- 

- 

- 

- 

- 

- 

7.1 

3 

- 

- 

- 

4.7 

- 

8.4 

7.1 

4 

- 

- 

- 

42.8 

- 

- 

- 

5 

- 

- 

66.1 

7 

8 

9 

- 

- 

- 

- 

- 

- 

5.3 

- 

- 

- 

- 

24.4 

- 

0.1 

10 

0.4 

0.4 

0.4 

0.4 

63.7 

0.4 

9.2 

11 

- 

- 

- 

- 

- 

- 

26.9 

12 

- 

- 

- 

- 

- 

- 

14.3 

13 

- 

- 

- 

- 

- 

- 

17,8 

TABLE  II. 4. 2 


III. 


IMPROVEMENTS  AND  MODIFICATIONS 


III.l  NORSAR  on-line  system  using  IBM  4331/4341  and  MODCOMP  Classic 
We  refer  to  the  detection  processor  operation  statistics  for  detailed 
performance  of  the  new  system. 

During  this  reporting  period  we  finished  array  monitor  control  (AMC) 
programs.  We  have  converted  the  programs  and  operational  routines  that 
were  used  for  the  old  SPS-system.  However,  the  system  has  now  capabilities 
for  monitoring  that  were  not  easily  accessible  using  the  360  system.  Using 
the  graphic  displays  we  may  now  do  extensive  visual  control  of  the  test 
period,  that  is,  we  may  for  instance  inspect  the  individual  broadband 
pulses  and  sine  wave  calibration  patterns  to  verify  the  quality  of  the 
instruments. 

The  IBM  2701  synchronous  controller  which  interfaces  IBM  4331  with  the 
MODCOMP  system  had  one  serious  breakdown.  It  was  only  due  to  our  own 
skilled  technician  that  the  system  was  brought  up  in  satisfactory  status 
again.  In  fact,  this  unit  has  been  functioning  with  even  fewer  bit¬ 
errors  during  data  transfer  after  this  breakdown.  Data  checks  (bit¬ 
errors)  on  this  unit  are  now  rarely  observed. 

We  have,  however,  had  occasional  data  losses  due  to  the  sharing  of  disk 
between  the  IBM  4331  and  4341.  Data  from  the  array  are  written  on  disk 
by  4331  and  are  read  from  the  same  disk  by  4341.  The  IBM  3370  disk  drives 
are  organized  in  so-called  strings.  On  one  string  we  had  disks  contain¬ 
ing  data  from  the  on-line  system  and  disks  containing  normal  user  area 
and  system  disks.  Since  the  4341  computer  is  significantly  faster  than 
the  4331  computer,  this  system  did  occasionally  block  the  slower  machine 
from  getting  access  to  the  disk.  The  follow-on  result  of  this  was  that 
the  disk  recording  on  4331  was  going  slower  than  the  data  was  coming 
in.  This  could  in  periods  cause  loss  of  up  to  a  minute  of  data.  We  solved 
this  problem  by  doing  two  things:  first,  we  reorganized  disk  strings  to 
separate  data  disks  from  system  disks;  second,  we  reorganized  programs 
on  4341  to  read  more  data  per  read  accesss.  Using  4341  and  3370  disks. 


a  single  read  operation  will  always  have  30  Msec  search  plus  data  trans¬ 
fer.  Reading  0.5  sec  of  data  will  hold  the  disk  for  33.4  Msec,  whereas 
reading  12.0  sec  of  data  will  hold  the  disk  for  63  Msec,  l.e.,  by  doubling 
the  time  we  keep  the  disk  busy,  we  may  get  24  times  more  data  to  work  with. 

The  result  of  these  operations  is  that  the  4341  uses  the  shared  disk 
significantly  fewer  times  per  time  unit.  By  applying  this  method  to  all 
programs  we  have  succeeded  in  completely  removing  this  form  for  loss 
of  data. 

During  the  second  half  of  this  reporting  period  we  also  started  using 
the  shared  disk  principle  for  the  programs,  that  is,  the  programs  run¬ 
ning  on  the  4331  computer  are  loaded  from  shared  disk.  We  therefore 
now  do  all  program  development/changes  on  the  4341  system,  and  a  new 
version  is  set  out  in  operation,  causing  less  than  2  minutes  drop  in 
data  recording. 

On  the  MODCOMP  system  we  have  experienced  that  the  7th  line  may  after 
a  few  days  of  continuous  operation  loose  synchronization.  We  may  cor¬ 
rect  this  by  commanding  the  MODCOMP  to  resynchronize  the  line  (by  in 
fact  doing  this  on  the  4341  system,  which  signals  4331  to  send  the 
command  to  MODCOMP),  or  restart  the  MODCOMP  system.  A  complete  restart 
of  the  MODCOMP  causes  less  than  one  minute  data  drop.  We  have  therefore 
preferred  to  restart  the  MODCOMP  system  up  to  twice  a  week  to  insure 
complete  resynchronization  of  all  lines  and  correct  timing.  This  effect 
is  a  result  of  the  SLEM  communication  being  critically  time  dependent. 

The  newer  (compared  to  SPS)  communication  interfaces  do  handle  some 
telephone  line  problems  better,  causing  less  number  of  transmission 
errors,  but  these  corrections  need  time  and  the  system  as  a  whole  may 
have  problems  in  keeping  all  seven  lines  exactly  synchronized.  We  have 
verified  this  by  dropping  subarray  beamforming  and  filtering  to  see 
whether  the  high  load  on  the  system  could  be  the  problem,  however,  this 
does  not  seem  to  be  a  problem.  In  conclusion,  we  do  not  see  any  data- 
loss  problem  here,  since  the  restart  or  resynchronization  twice  a  week 
will  secure  the  overall  performance  of  the  system. 


The  implementation  of  an  IBM  Series  1  front-end  processor  has  given 
the  opportunity  to  dial  the  NORSAR  computer  center  from  a  remote 
location.  Using  a  cheap,  standard  ASCII-type  terminal,  we  may  dial 
the  center  and  use  the  facilities  of  the  4341  system.  This  has  been 
used  by  DP  operators  to  check  the  status  of  the  system.  In  fact,  we 
have  verified  that  this  may  also  be  done  from  Washington,  D.C. ,  or 
any  other  place,  using  special  300  baud  modems. 

Improvements  of  the  system  now  would  be  to  exchange  the  2701  communica¬ 
tion  controller  with  a  direct  IBM  channel  interface  device  on  the  MODCOMP. 
The  main  reason  is  to  avoid  having  a  device  that  is  20  years  old,  and 
which  in  the  near  future  no  longer  will  be  supported  by  IBM.  (However, 
NORSAR  is  glad  to  have  a  technician  who  is  able  to  repair  components  on 
such  devices.)  Another  gain  we  would  experience  by  doing  so  is  that 
the  data  transfer  time  between  IBM  and  MODCOMP  will  be  reduced  by  a 
factor  of  5,  which  will  free  some  of  the  load  in  both  ends. 

The  next  improvement  step  would  be  to  build  up  backup  systems.  Today 
we  may  run  the  IBM  part  of  the  system  on  either  4331  or  4341,  but  there 
is  no  backup  for  MODCOMP.  There  are  several  wayB  to  do  this:  One  is 
of  course  to  have  an  identical  MODCOMP  system  as  backup.  Another  is 
to  move  subarray  beamforming  and  filtering  over  to  the  IBM  side  (if 
capacity).  This  reduces  the  requirements  for  the  MODCOMP  backup  machine, 
and  increases  the  possibilities  for  testing  new  DP  processors  on  the 
IBM  side.  This  latter  way  reduces  the  MODCOMP  system  to  a  pure  communi¬ 
cation  path  to  the  subarrays.  This  is  important  since  the  next  step  in 
future  Improvements  of  the  NORSAR  array  would  be  on  the  instrument  side, 
that  is,  the  SLEM  itself  might  be  considered  replaced. 


J.  Fyen 
T.  Hoff 
R.  Paulsen 
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Appendix  to  III.l  -  Program  examples,  corrections  to  NORSAR  Scientific 
_ Report  No.  1-82/83 _ 


The  following  are  parts  of  a  FORTRAN  G1  program  that  can  extract  data 
from  NORSAR  high-rate  data  block. 


L 


C 


C 


LOGICAL* 1  NRDATA (13721 

LOGICAL*  1  NRIO  (?)  .NRTIHE  (A)  , NRGEM 2  )  ,  AUT  S  i  I  ?  )  •  AUTS  E  1(21) 
LOGICAL*!  HANSAI2) .PARSE  I (21)  ,RULTI l 7 ) , PAO .LOGR II 262 ) .L0GR2 ( 242  I 
LOGICAL* 1  L0GR3I262)  . LCGR4 1 262 ) , L 0GR5 I  262 1 . LOGR I  1310) 

EOL'I  VALENCE  (NROA  TA .NRI  0  )  .(NRDATA  (3  )  .  NRT  IRE  1  .  ( NRCAT  A  (  7  )  .NRCEN  I , 

•  (NR0ATAI9) • A  UT  SA ) • ( NROA T  A (11 ) , AUT  SE1 ) .INRDATAI32). MANS A ) . 

•  (NRDATA (34  I .PANSEI I , (NROATA(SS) .MULTI ) , (NROATAI 62 ) .PAD ) • 

•  (NRDATA  (63)  *  L0GR1-.LCGR  )  .  ( NRO  A  T  A  (  3  25 )  ,  LOG  R2  )  .  ( NRC4T  At  5(7).  LOGR  3). 

•  (NRDATA l 849)  . LCGR4 ) . ( NRO AT  A ( 1 1 1 1 )  .LCGR9 ) 

LOGICAL* 1  LOGREC (262  I .STATUS (10) .SOBL (18 .2,7 ) 

IN  TEGER*  2  SOB  (  9  >2 ,7  ) 

E0L1 VALENCE (LOGREC. STATUS  I . (LOGREC ( 1 1  I  ,S CEL . SCE ) 


ONLCC740 

CNICC1!C 

ONLOOTtO 

ONL00770 

ONLCOTIO 

CNL00C10 

ONL00C20 

ONLOO03O 

CNL00E4C 

GNLCCEiC 

ONLOC0EO 

ONLOO07O 

ONLCC08C 

CNLCC09O 


C 

R E AO  THE  EQUIVALENT  OF  0.5  SECCNCS  CATA  (1372  0VTES  NRCATA) 

C 

C  PICK  UP  SINGLE  SENSORS,  SR  PROP  NORSAR  ARRAY 

C 

C  N2C4C  «  1  FOR  10  HZ 
C  N2C40  •  2  FOR  20  H2 

00  3C00  ILOG-1.5 

C 

ILDIJP  * ( I  LOG-1 ) *262 

C 

00  2100  1*1, 242 

2 1 C  0  LOGREC  ( I  )  •  LOGR ( I* I LD I SF) 

C 

00  2430  120*1, N2J40 

C 


00  2400  I  SUB*  1,  7 

C 

ICWOOW  •  SOS (1,120.1  SUB  I 

c 

C  1 F (  ICWOOW  . NE.  32  )  GO  TO  ERROR 

C 

00  2320  ISNS*  1,  6 

C 

ISANP  •  I SAHP+1 

SAHP ( I  SAMP )  •  SOB (3*ISNS« 120, ISU6  )/4 


C 

C 

2320  CONTINUE 

C 

c 

24C0  CONTINUE 

C 

c 

j;:c  CONTINUE 
c 


ONL02000 
0NLC2C1C 
ONL02C20 
ONL02020 
ONLO2030 
ONL02C40 
CNL02C5C 
ONL02C6C 
ONLO2O0O 
ON  L  0  2  C  9  G 
CNL0210C 
ONLO2110 
ONL02120 
ONL0213C 
CN102K0 
ONL02170 
ONLO210O 
0NL021SC 
CN1022CC 
□NL0221C 
ONL02230 
ONL  0224C 
CNL0225C 
ONL0224C 
0NLO2270 
ONLO220O 
0NL02300 
0NL02310 
ONL02320 
ONL0234C 
ONIO230C 
ONL023A0 
ONLO240O 
ONL024SC 
CNLOiioe 


C  CK102C4C 
C  LOW  RATE  SENSOR  DATA  EXTRACTION  ONLQ204O 
C  SEE  NORSAR  SC1ENTIF1  REPORT  NO.  1-B2/S3  ONL02440 
C  ONLOICSO 


III. 2  Description  of  NORSAR  recording  system 

Over  Che  last  few  years,  the  number  of  independent  recording  system  used 
at  NORSAR  has  increased  steadingly,  and  with  that  the  need  for  easy  and 
standardized  retrieval  of  information  about  response  characteristics  and 
quantification  levels.  This  has  now  been  done  in  the  form  of  a  subroutine 
from  which  the  desired  information  is  returned  following  the  entering  of 
system  number  and  frequency.  The  call  is  as  follows: 

CALL  SPSYST  (J,  FREQ,  RESP,  GAIN,  SENS,  CODE,  NSYS,  IERR) 

where  input  is 

J  *  System  number 

FREQ  -  Frequency  (Hz)  for  which  response  values  are  wanted 
and  output  is 

RESP  -  Response  value  relative  to  1  Hz 
GAIN  -  Response  in  NM/QU  at  1  Hz  for  digital  systems 
SENS  »  Absolute  magnification  at  1  Hz  for  analog  systems 
CODE  ■  System  name  (REAL  *  8) 

NSYS  «  Number  of  system  responses  on  file 
IERR  -  Error  code,  0  -  ok,  1  otherwise 

The  data  base  consists  presently  of  19  systems  as  described  in  Table  III. 2.1, 
where  much  relevant  Information  is  given,  except  for  specifications  about 
amplifiers  and  AD-convertors.  Only  two  of  the  systems  (5-6)  are  analog, 
based  on  MEQ-800  recorders.  It  should  be  emphasized  that  phase  information 
is  not  given  for  any  of  the  systems. 
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Table  III. 2.1  Description  of  the  19  recording  systems  used  at  NORSAR  over  the  last  few  years, 
with  recording  unit,  amplification,  seismometer  type,  anti-aliasing  filter, 
quantification  at  1  Hz,  and  comments. 
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The  standard  NORSAR  system  (1,2,15)  is  gain-ranged  in  the  following  way: 


antum  levels 

Steps 

NM/QU 

1-127 

1 

0.0427 

128-511 

4 

0.1708 

512-2047 

16 

0.6832 

2048-8191 

64 

2.7328 

The  other  system  with  gain-ranging  is  the  Kinemetrics  PDR-2,  where  there 
is  a  12  bit  mantissa  (±2047  QU)  and  a  3  bit  Log  (base  2)  gain  value  with 
7  steps  from  2®  (0  dB)  for  the  strongest  signals  to  T?  »  128  (42  dB)  for 
the  weakest  signals.  The  gain  switching  logic  is  such  that  the  system  is 
slow  to  amplify  but  fast  to  attenuate. 

The  relative  response  values  for  different  frequencies  for  all  the  systems 
are  given  in  Table  111.2.2,  with  3  values  per  octave  in  frequency.  These 
are  the  values  that  are  on  file  (in  SPSYST),  and  values  for  other  fre¬ 
quencies  are  obtained  by  linear  interpolation  between  the  two  closest 
points. 

The  relative  response  functions  are  plotted  in  Fig.  III. 2.1  for  all  of  the 
systems  and  in  Fig.  III. 2. 2  for  the  systems  that  have  been  used  more  fre¬ 
quently  lately.  The  curves  are  mostly  based  on  theoretical  values,  although 
most  of  them  are  also  checked  by  independent  calibration.  It  should  be  em¬ 
phasized  that  the  wide  frequency  range  has  been  chosen  for  reasons  of  con¬ 
venience,  values  on  the  flanks  have  been  extrapolated  and  are  therefore 
not  necessarily  realistic. 

Weathejr  £tatio£s_(wirui_direc£i£n_and_8j>e£d_) 

Two  weather  stations  have  recently  been  installed  at  NORSAR,  one  at  NORESS 
(06C)  and  one  at  subarray  01A. 


6  3114 
M31SAS 


Table  III.2.2  Systems  response  values  for  the  19  systems  described  In  Table  III. 2.1.  The 

values  are  logarithmic ,  and  correspond  to  the  curves  plotted  in  Figs.  III.2.1-! 


RELATIVE  DISPLACEMENT 


SYSTEM  RESPONSE  FUNCTIONS 


FREQUENCY  (HZ) 


Fig.  III.  2.1  System  response  for  all  19  systems  di 
section.  The  system  numbers  refer  to 
Tables  III.  2. 1-2. 
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SYSTEM  RESPONSE  FUNCTIONS 


1. 

2. 

9. 

12. 

13. 

14. 

15. 

16. 

17. 

18. 
19. 
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Fig.  III. 2.2  System  response  for  the  systems  more  recently  used  in  the 

analysis  of  NORSAR  and  NORESS  data.  The  system  numbers  refer 
to  the  numbers  In  Tables  III. 2. 1-2. 
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The  wind  data  from  NORESS  are  transmitted  in  analog  form  to  the  Modcomp 
at  NDPC  (system  13  in  Table  III. 2.1),  with  wind  direction  in  channel  7 
and  wind  speed  in  channel  8.  At  subarray  01A,  equivalent  data  are  recorded 
at  channels  01A05  and  01A03,  respectively.  The  channels  are  calibrated 
a^  follows : 


Channel 

Data 

5V  level 

mV/QU 

Quantification 

NORESS  7 

Direction 

540° 

3.40 

0.367  deg/QU 

NORESS  8 

Speed 

50  m/s 

3.40 

0.034  m/s/QU 

01A05 

Direction 

540° 

0.61 

0.066  deg/QU 

01A03 

Speed 

50  m/s 

0.61 

0.0061  m/s/QU 

Note  here  that  the  wind  speed  Increases  linearly  from  0  to  50  m/s 
corresponding  to  0  to  5  volts  in  output,  while  the  same  output  range 
for  the  wind  direction  corresponds  to  a  linear  clockwise  increase  from 
0°  (north)  to  540°  (south).  A  further  increase  will  cause  the  output 
to  drop  from  5  volts  ■  540°  to  5/3  volts  «  180°,  while  a  counter¬ 
clockwise  change  around  0°  will  cause  the  output  to  increase  from  0  volts 
-  0°  to  10/3  volts  »  360°.  This  is  done,  obviously  enough,  in  order  to 
reduce  the  number  of  sudden  changes  in  the  output  level  as  the  wind 
direction  fluctuates  around  north  or  south. 


H.  Bungum 


IV.  FIELD  INSTRUMENTATION  AND  MAINTENANCE  ACTIVITIES 
Improvements  and  modifications 


Reference  is  made  to  Table  IV. 1  and  IV. 2  indicating  the  status  of  the  SP 
instruments  (original  array),  the  modified  NORESS  array,  and  the  expanded 
02B  subarray  (telemetry  stations),  which  are  all  unchanged  compared  to  status 
in  the  previous  report. 

The  expanded  02B  subarray  was  fully  operational  5  November  1982. 

Fig.  IV. 1  shows  the  expanded  02B  array,  and  Table  IV. 3  gives  the  coordinates 
for  the  stations. 

The  06C  NORESS  array  is  shown  in  Fig.  IV. 2. 

Preparatory  work  in  connection  with  the  planned  expansion  of  the  NORESS  array 
continued  at  the  maintenance  center  (NMC)  in  March. 

In  connection  with  a  noise  study  program,  recording  of  noise  data  started 
early  in  January  and  continued  throughout  the  reporting  period,  reference 
Table  IV. 4. 
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Station 

No. 

Latitude 

(N) 

Longitude 

(E) 

1 

61.068344 

11.156468 

2 

61.101196 

11.161124 

3 

61.091309 

11.166824 

4 

61.107315 

11.174083 

5 

61.084824 

11.174442 

6 

61.049728 

11.158080 

TABLE  IV. 3 

Preliminary  geographical  coordinates 
for  the  new  02B  stations. 


Fig.  IV. 1  The  six  new  stations  in  the  expanded  02B  array 


Activities  in  the  field  and  at  the  NORSAR  Maintenance  Center 


This  section  outlines  in  brief  the  activities  in  the  field  and  at  the 
NORSAR  Maintenance  Center  (Table  IV. 4).  Table  IV. 5  indicates  corrective/ 
preventive  maintenance  in  the  field. 


SA/Area 

Task 

Date 

01A 

Installation  of  a  wind  station 

(Week  ..  Oct) 

01A 

SP  seism,  installed  in  borehole 

(7,29  Oct) 

06C,  NORESS 

Inst,  of  two  sets  of  3-comp 

SP  seism.  (SS-1)  6  a  wind  station 

(1,4, 5, 6  Oct) 

02B,  telemetry 

sta. 

Inst,  of  equipment  at  the  new 
stations  continued 

(13,14,18,19, 
20  Oct) 

SNSN,  Evje 

Equipment  repaired 

(22,23  Oct) 

02B,  telemetry 

sta. 

Inst,  of  equipment  finished  and 
fully  operational 

(5  Nov) 

06 C,  NORESS 

Attempts  to  locate  new  sites 

(10,11  Nov) 

02B 

Subarray  visited  due  to  power  outage 

(20  Dec) 

02B,  telemetry 

sta. 

Snow/ice  problems  reducing  charging 
effect  of  solar  cells 

06C,  NORESS 

Array  expansion 

(1,15,16  Dec) 

NORSAR 

Meeting  at  NORSAR,  Sandia  repres. 

(6  Dec) 

06C  area 

P.W.  Larsen  guided  Sandia  and 

NORSAR  repres.  in  connection  with 
new  array 

(8  Dec) 

03C  area 

Assistance  to  Rena  Power  Co.  in  con¬ 
nection  with  power  outage 

(3,4  Jan) 

03C  area 

Noise  data  recording 

(11  Jan) 

02B  area 

(12  Jan) 

02C  area 

— H  — 

(19  Jan) 

04C  area 

_ll_ 

(25  Jan) 

02B,  telemetry 

sta. 

Sta.  1,4  -  battery  replaced 

(14,20  Jan) 

02B  area 

Noise  data  recording 

(1,25  Feb) 

06C,  CTV 

T.O.D.  generator  hooked  up  for 
M0DC0MP  timing  comparison 

(7  Feb) 

06C  area 

Noise  data  recording 

(8,22  Feb) 

06 C  area 

Weather  station  checked 

(9  Feb) 

03C  area 

Localization  of  cable  traces 

(16,18  Feb) 

01 A  area 

Noise  data  recording 

(24  Feb) 

02B  area 

(1  Mar) 

02B,  telemetry 

sta. 

Noise  problem  investigation  cont. 

(2  Mar) 

01A  area 

Noise  data  recording 

(3,28  Mar) 

NMC 

Activities  mainly  in  connection 
with  preparation  of  the  planned 
expansion  of  the  NORESS  array 

(4-23  Mar) 

06C  area 

Noise  data  recording 

(24,25  Mar) 

NDPC 

(  ) 

TABLE  IV. 4 


Period  1  Oct  82  -  31  Mar  83 


Unit 

SLEM  Equipm. 

SP  Equipm. 

LP 

Equipm. 

Rep.  Adj.  Repl. 

Rep.  Ad.i.  Repl. 

Rep. 

Adj.  Repl. 

Seism.  LP 

FP 

HP 

ROD,  FP 

ROD,  MP 
MP-bridge 

1 

1 

3 

5 

1 

Ampl .  Ithaco 

LP 

1 

Ampl .  RA-5 

SP 

1 

Seism.  SP 

Vert.  (SS-1) 
NORESS 

1 

SLEM 

Mux/ADC 

1 

TABLE  IV.  5 

Corrective/preventive  maintenance  work  in  the  period. 
1  Oct  82  -  31  Mar  83 


Array  statue 

As  of  31  March  1983  the  following  channels  deviated  from  tolerances: 
01B01 ;  02B03.05;  02006,08,09;  03C04;  04C01;  06C09. 


01 A 

01 

8  Hz 

filter 

01A 

02 

,  60  m  borehole 

01A 

03 

Wind 

speed  measurements 

01A 

04 

Attenuated  40  dB 

01A 

05 

Wind 

direction  measurements 

0. A.  Hansen 
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ABBREVIATIONS 


CTV 

- 

Central  Terminal  Vault 

EW 

- 

East-West 

FP 

- 

Free  period 

LP 

- 

Long  period 

HP 

- 

Hass  position 

MUX 

- 

Multiplexer 

NDPC 

- 

NORSAR  Data  Processing  Center 

NORESS 

- 

NORSAR  Experimental  Small-Aperture  Subarray 

NS 

- 

North-South 

RCD 

- 

Remote  centering  device 

SLEM 

- 

Seismic  short  and  long  period  electronics  module 

SA 

- 

Subarray 

aP 

- 

Short  period 
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VI.  SUMMARY  OF  TECHNICAL  REPORTS /PAPERS  PREPARED 
VI. 1  The  absorption  band  effect  on  long-period  body  waves 
It  is  generally  appreciated  that  the  effect  of  anelastic  absorption  on 
body  wave  amplitudes  is  most  pronounced  in  the  short-period  band.  However, 
in  the  presence  of  an  absorption  band  in  the  mantle,  the  accompanying 
velocity  dispersion  is  best  observed  in  long-period  body  waves.  This  is 
because  (1)  dispersion  represents  an  integral  form  of  Q”1  in  the  band  (in 
contrast  to  attenuation  which  represents  a  'point  property'),  and  (2)  the 
mantle  absorption  band  extends  throughout  the  long-period  range,  but  not 
throughout  the  short-period  range  of  body  waves.  The  last  statement  is  an 
inference  drawn  by  several  authors  in  recent  years,  and  it  is  confirmed 
in  the  present  study.  In  this  work  the  possibility  is  explored  to  use 
both  the  dispersion  and  the  amplitude  information  in  a  procedure  to 
separate  absorption  and  source  effects,  and  the  procedure  is  applied  to 
SRO/ASRO  records  of  P  waves  from  moderately  sized  deep  focus  events 
between  30  and  90°.  The  observations  were  compared  to  a  synthetic 
purely  elastic  response,  and  the  observed  short-period  onset  was  used  to 
adjust  the  start  time  of  the  synthetic.  The  last  procedure  eliminates 
the  effect  of  station  residuals.  The  implicit  assumption  here  is  that 
the  visible  short-period  onset  is  unaffected  by  the  absorption;  this  can 
be  justified  later  by  the  results  of  the  analysis.  The  data  are  shown 
in  Fig.  VI. 1.1;  other  information  can  be  found  in  Doornbos  (1983). 

In  all  examples  the  long-period  waveform  is  essentially  delayed  with 
respect  to  the  synthetic,  and  the  peak  amplitude  delay  is  used  as  an 
observational  parameter.  The  delay  must  be  explained  as  the  combined 
effect  of  the  source  pulse  and  the  anelastic  dispersion.  Since  dispersion 
is  relatively  Insensitive  to  details  of  an  absorption  model,  we  choose  a 
convenient  model:  the  relaxation  model  which  is  uniform  in  Jtnr,  within 
certain  bounds  Tj,  13  (Liu,  Kanamori  and  Anderson,  1976): 

w(t2~Ci) 

Q-l(u)  ~  2  tg"1  { - j  (1) 

1-Hd2tjT2 
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Our  measure  of  phase  velocity  dispersion  Is: 

D(u>)  «  c„c(u>)-l  *  Q~1(u)  (2) 

where  Cm  and  c(<i>)  are  phase  velocities  at  infinite  frequency  and  at  ui, 
and  the  overbar  denotes  the  Hilbert  transform*  For  the  model  (1): 

1+w2T2 

D(<*>)  ~  ln(T2/T-.)-i  4n( - — )  (3) 

1+<i)2t^ 

This  model  is  convenient,  since  the  resulting  damping  and  dispersion  become 
explicitly  Independent  of  the  long-period  cut-off  T2»  for  WT2»1.  If 
also  T2»t^: 

Q-1(u>)  =  2 h  ^tg^d/csxj),  D(w)  *  */2  in(l+l/u>2tj) 

There  is  good  evidence  that  UT2»1  at  body  wave  frequencies  (Anderson 
and  Given,  1982).  The  second  condition  T2»Tj  can  be  justified  a  posteriori. 
It  sets  the  normalizing  constant  in  equation  (4)  to  x.  In  the  present  treat¬ 
ment  we  have  further  specified  the  model  by  assuming  that  is  given  by 
PREM  (which  incorporates  a  long-period  Q  model).  There  is  then  one  undeter¬ 
mined  parameter  t}. 

In  treating  the  source  effect  two  approaches  can  be  taken.  The  first  is  to 
ignore  the  effect,  i.e.,  treat  the  source  pulse  as  a  delta  function.  It 
results  in  a  maximum  estimate  for  the  cut-off  frequency  of  the  absorption 
band,  i.e.,  an  absolute  minimum  for  Tj.  In  the  second  approach  we  adopt  a 
model  for  the  source  pulse  and  estimate  the  source  and  absorption  parameters 
simultaneously.  It  is  done  by  computing  the  RMS  amplitude  error,  and  the 
time  residue,  as  a  function  of  trial  cut-off  relaxation  time  tj. 

The  results  of  this  experiment  are  shown  in  Fig.  VI. 1.2.  The  time  residue 
and  RMS  error  curve  are  consistent,  and  zero  time  residue  corresponds  to  an 
absorption  bandwidth  near  0.2  Hz,  i.e.,  near  0.9  s.  In  the  presence  of 
noise  and  in  the  context  of  a  single  absorption  band  this  constraint  can  be 
somewhat  relaxed,  and  cut-off  frequencies  up  to  1  Hz  are  admitted  by  the 


data.  The  result  represents  the  integrated  effect  of  Q“1  along  the  wave 
path  but  it  can  be  directly  applied  to  correct  body  waves  synthesized  in 
purely  elastic,  or  frequency  independent  Q,  models.  Since  the  present 
result  is  based  on  teleseismic  records  from  deep  events,  it  should  be 
more  representative  of  absorption  in  the  lower  mantle. 


D.J.  Doornbos 
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Fig.  VI. 1.2  Averaged  time  residue  and  RMS  amplitude  error  In  fitting  a 

source  and  absorption  model  to  the  P  wave  data  of  Fig.  VI. 1.1 
as  a  function  of  the  high-frequency  cut-off  of  the  absorption 
band.  The  closed  dots  represent  a  solution  with  finite  source 
effect.  The  open  dots  represent  a  solution  where  the  source 
effect  Is  Ignored. 


VI. 2  The  present  seismic  evidence  for  a  boundary  layer  at  the  base  of 
the  mantle 

The  velocity  structure  at  the  base  of  the  mantle  has  been  the  subject  of 
a  long-standing  controversy.  Velocity  models  proposed  in  the  past  encompass 
a  range  from  positive  to  negative  gradients;  the  three  models  shown  in 
Fig.  VI. 2.1  are  representative  of  different  alternatives.  This  figure 
also  shows  two  different  interpretations  of  combined  P  and  S  velocity 
profiles,  as  a  density  effect  or  as  a  temperature  effect.  Note  that  models 
with  velocities  which  are  anomalous  in  the  sense  that  the  gradients  with 
depth  reverse  and  become  negative,  have  recently  been  associated  with  a 
thermal  boundary  layer.  In  1974,  Cleary  reviewed  the  available  seismic 
evidence  and  concluded  that  a  velocity  reversal  was  likely.  However,  the 
data  were  interpreted  mainly  in  a  ray  geometrical  context  and  since  dif¬ 
fraction  effects  can  be  significant,  the  evidence  has  been  subject  to 
criticism.  Doornbos  and  tiondt  (1979)  applied  diffraction  theory  to  an 
extensive  set  of  long-period  P  and  SH  data  and  still  inferred  a  velocity 
reversal,  although  the  thickness  of  the  low-velocity  zone  was  found  to 
be  much  smaller  than  the  200  km  conventionally  assign  to  the  D"  layer. 

More  recently,  several  authors  find  no  evidence  for  a  low-velocity  zone 
(see,  e.g.,  Ruff  and  lielmberger,  1982).  In  a  recent  contribution  we  have 
made  use  of  new  data  from  the  SRO  network  and  from  the  NORSAR  array, 
in  an  attempt  to  reconcile  the  available  evidence. 

A  typical  example  of  long-period  diffracted  waves  observed  at  the  SRO 
network,  is  shown  in  Fig.  VI. 2. 2.  It  is  clear  that  the  amplitude  decay 
of  SV  in  the  core  shadow  zone  is  too  strong  to  yield  useful  data,  in 
agreement  with  a  theoretical  prediction  for  a  range  of  possible  models. 
Thus,  only  P  and  SH  have  been  used  here.  Some  striking  examples  of  short- 
period  diffracted  P  at  the  NORSAR  array  are  given  in  Fig.  VI. 2. 3.  PKP 
from  the  same  events  are  given  for  comparison,  and  the  data  to  be  used 
are  based  on  an  amplitude  ratio  of  P  and  PKP  at  the  periods  1  and  2 
seconds.  Altogether  we  have  obtained  long-period  diffraction  data  from 
16  events,  long-period  P  from  25  events,  and  short-period  P  from  10 


events 
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In  Figs*  VI. 2. 4  and  VI. 2. 5  the  data  are  compared  to  calculated  decay 
spectra  and  dT/dA,  using  the  models  of  Fig.  VI. 2.1.  In  the  calculations, 
the  extended  WKBJ  method  (Langer's  approximation)  was  applied  to  a 
piecewise  smooth  layered  model.  In  this  way  the  effect  of  both  a  velo¬ 
city  gradient  and  changes  in  the  gradient  are  evaluated  (Doornbos,  1981). 
The  dT/dA  curves  of  Fig.  VI. 2. 4  reflect  the  fact  that  diffracted  waves 
are  usually  dispersed.  In  analogy  to  phase  and  group  velocity,  we  have 
introduced  the  phase  and  group  slowness,  p(ui)  and  g(u),  respectively. 

Their  difference  is 

g(w)  -  p(w)  ■  is  dp/du) 

and  since  this  dispersion  depends  on  the  velocity  gradient  and  its  changes, 
it  can  be  used  as  an  additional  diagnostic. 

The  observed  low  decay  of  long-period  diffracted  SH,  the  dT/dA  values,  and 
the  apparent  absence  of  dispersion,  are  all  compatible  with  a  relatively 
thin  (<  100  1cm)  low-velocity  zone  for  S  at  the  base  of  the  mantle.  The 
long-period  SRO  data  lack  the  resolving  power  to  settle  the  question  for 
P.  The  results  for  short-period  P  as  summarized  in  Fig.  VI. 2. 5  strongly 
suggest  that  a  comparable  low-velocity  zone  for  P  exists,  at  least  in  the 
region  sampled  by  the  NORSAR  data  (beneath  Central  Asia).  It  is  not  in¬ 
consistent  with  a  global  average  of  short-period  dT/dA  and  amplitude  data 
near  the  core  shadow  boundary  (also  summarized  in  Fig.  VI. 2. 5),  but  it 
should  be  mentioned  that  at  least  one  recent  study  of  short-period  ampli¬ 
tudes  gave  results  conflicting  with  the  global  average  (Ruff  and  Helmberger 
1982).  If  this  is  to  be  'explained'  by  lateral  variations,  then  it  should 
be  added  that  the  presence  of  lateral  variations  itself  is  not  in  conflict 
with  boundary  layer  models.  On  balance,  current  thermal  models  appear  to 
be  consistent  with  (but  not  required  by)  the  present  seismic  evidence. 


D.J.  Doornbos 
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Fig.  VI. 2.1  Models  of  P  and  S  velocity  Vp,vs,  density  p  snd  temperature  T 
In  the  lower  mantle.  The  super diabetic  temperature  profiles 
of  PREM  and  LI  correspond  to  a  homogeneous  density  (PEM 
profile),  and  the  inhomogeneous  density  profiles  correspond 
to  the  adiabatic  temperature  of  PEM.  Absolute  temperatures 
have  been  arbitrarily  fixed. 
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Fig.  VI. 2. 2  Vertical,  radial  and  transverse  component  SRO  records  of  P,  SV,  and  SH  near 
a  great  circle  path  around  the  core.  Time  length  Is  2  minutes. 
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Fig.  VI. 2. 4  Decay  spectra  and  dT/dA  of  P  and  SH.  The  decay  y  measures 
logarithmic  amplitude  decay  per  degree  of  eplcentral  dis¬ 
tance  and  corrected  for  geometrical  spreading.  Group-phase 

slowness  is  a  measure  of  dispersion.  Frequency  point  1 . 8 

correspond  to  0. 015625, .... ,0.125  Hz  in  steps  of  0.015625. 

Theoretical  curves  for  models _ :  PEM; _ :  PREM; 

_  :  LI.  Averaged  observational  data  •  :  this  study; 

□  and  V  :  other  studies. 
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Fig.  VI. 2. 5  dT/dA  and  relative  amplitudes  of  short-period  P.  Theoretical 

curves  for  node Is  PEM,  PREM  and  LI.  Amplitude  bands  correspond 
to  period  range  1-2  seconds.  Additional  aaplltude  decay  due 
to  anelaatlc  damping  Indicated  for  aodel  LI  with  Qa  »  400  at 
base  of  aantle.  taplitude  data  between  115-135°  are  froa  NORSAR 
(this  study).  Aaplltude  and  dT/dA  data  between  80-105°  are  froa 
other  studies. 
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VI. 3  The  effect  of  aliasing  on  the  NORSAR  detector  performance 
The  standard  NORSAR  data  are  recorded  with  a  sampling  rate  of  20  Hz, 
but  with  an  anti-aliasing  filter  at  5  Hz  (4.75  Hz  to  be  precise). 

The  reason  for  this  is  that  the  on-line  Detection  Processing  (DP) 
has  always  been  done  at  a  10  Hz  sampling  rate,  and  the  filter  has  been 
kept  low  enough  to  avoid  aliasing  problems  there. 

With  the  increased  interest  for  regional  phases  and  higher  frequencies 
the  question  has  been  raised  if  the  filter  could  be  changed  from  5  Hz 
to  S  Hz,  and  we  have  therefore  now  investigated  the  effect  of  this  on 
the  on-line  detector  system.  This  has  been  done  on  the  basis  of  data 
that  have  been  recorded  experimentally  with  8  Hz  filters  at  some 
individual  channels.  Host  of  the  events  analyzed,  shown  in  Table  VI. 3.1, 
are  taken  from  a  time  period  in  1978/79  when  we  had  5  subarray  center 
seismometers  equipped  with  8  Hz  filters  (01A00,  02B00,  02C00,  04C00, 
06C00),  and  the  rest  of  the  data  are  from  1982/83,  with  only  2  '8  Hz 
channels'  in  operation. 

We  have  selected  events  (see  Table  VI. 3.1)  from  regions  with  particularly 
high-frequency  signals,  which  are  the  signals  that  presumably  should  be 
most  affected  by  aliasing  problems  as  a  result  of  combining  8  Hz  filters 
with  10  Hz  processing.  The  analysis  procedure  has  been  as  follows: 

1.  The  17  events  in  Table  VI. 3.1  were  selected  and  the  data  put  on 
disk  files. 

2.  For  each  event,  the  data  were  then  plotted  as  shown  on  the  example 
in  Fig.  VI. 3.1,  with  2  sampling  rates  (20  and  10  Hz)  and  2  filters 
(1.2-3. 2  and  2. 0-4.0  Hz).  These  filter  bands  are  identical  to  those 
presently  used  in  NORSAR  DP.  For  comparison,  the  unfiltered  trace 
was  also  plotted. 


3.  For  each  of  the  4  filtered  traces,  STA/LTA  was  computed  in  a  way 
as  close  as  possible  to  what  is  done  in  DP,  with  options  to  plot 
both  ST A,  LTA  and  STA/LTA  traces.  The  latter  ones  are  shown  in 
Fig.  VI. 3. 2  for  the  data  plotted  in  Fig.  VI. 3.1.  It  is  seen  from 
that  figure  that  we  do  lose  some  SNR  as  a  result  of  the  decimating 
process,  for  the  1 .2—3.2  Hz  the  SNR  drops  from  3.97  to  3.83  and  for 
the  and  for  the  2. 0-4.0  Hz  filter  from  7.27  to  6.40. 

The  results  for  all  of  the  events  analyzed  are  shown  in  Figs.  VI. 3. 3-4, 
plotted  individually  for  teleseismic  and  regional  events,  and  for  each 
of  the  2  filters  for  each  event  group.  We  can  draw  the  following  con¬ 
clusions  from  the  figures: 

1.  For  teleseismic  events  processed  at  1.2-3. 2  Hz,  a  change  from  5  to 
8  Hz  anti-aliasing  filters  would  not  affect  the  performance  of  the 
10  Hz  DP. 


2.  For  teleseismic  events  processed  at  2.0-4. 0  Hz  there  is  a  loss  in 
SNR  typically  around  1-2  dB.  For  the  events  analyzed  here  (and 
generally  for  most  Eurasian  events  recorded  at  NOKSAR)  this  filter 
band  has  a  significantly  better  SNR,  and  a  filter  change  as  dis¬ 
cussed  can  therefore  not  be  implemented  without  a  clear  negative 
effect  on  the  DP  performance. 

3.  For  regional  and  local  events  we  see  t.jai  Fig.  VI. 3. 4  that  the  scatter 
in  the  SNR  values  is  much  larger  than  for  the  teleseismic  events, 

but  with  no  systematic  difference  in  the  SNR  range  critical  for 
detections.  A  filter  change  should  therefore  not  be  a  problem  for 
these  events. 

The  above  analysis  has  been  done  on  single  channels,  and  it  should  be 
expected  that  the  high-frequency  beamforming  loss  should  decrease  a 
little  the  adverse  effects  of  a  filter  change  to  8  Hz.  Our  conclusion 
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is  nevertheless  that  an  overall  change  from  5  to  8  Hz  anti-aliasing 
filters  should  not  be  implemented  unless  combined  with  a  change  from 
10  to  20  Hz  processing  in  DP. 


0.  Pettersen 
H.  Bung urn 
F.  Ring dal 


No. 

Date 

Region 

Arr.  Time 

(deg) 

mb _ 

NCH 

Comments 

1 

77/12/26 

E.  Kazakh 

04.10.17 

38 

4.9 

5 

Expl . 

2 

78/05/29 

E.  Kazakh 

05.04.17 

38 

4.7 

4 

Expl. 

3 

78/09/20 

E.  Kazakh 

05.10.17 

38 

4.2 

3 

Expl . 

4 

78/12/20 

E.  Kazakh 

04.40.17 

38 

4.7 

4 

Expl . 

5 

79/09/06 

C.  Siberia 

18.07.17 

38 

4.9 

3 

Expl. 

6 

79/09/17 

S.  Norway 

15.22.18 

3 

2.8 

5 

Expl . 

7 

79/09/18 

S.  Norway 

17.56.09 

4 

2.6 

4 

Expl . 

8 

79/09/23 

N.  Norway 

13.54.29 

6 

3.1 

5 

9 

79/10/02 

S.  Norway 

14.19.15 

3 

2.4 

4 

Expl . 

10 

79/10/07 

C.  Siberia 

21.08.06 

44 

4.9 

4 

Expl. 

11 

82/12/25 

E.  Kazakh 

04.30.27 

39 

4.7 

2 

Expl . 

12 

83/01/05 

Honshu,  Japan 

07.29.44 

71 

5.1 

1 

13 

83/01/07 

Honshu,  Japan 

18.30.32 

74 

5.1 

1 

14 

83/02/13 

N.W.  Kashmir 

02.00.55 

47 

5.1 

2 

15 

83/02/18 

Persian  Gulf 

07.48.25 

47 

5.0 

1 

16 

83/02/20 

Hokkaido,  Japan 

01.32.44 

70 

5.1 

1 

17 

83/02/28 

Aleutian  Isl. 

05.44.23 

70 

5.2 

1 

Table  VI. 3.1  List  of  events  used  In  analysis  of  8  Hz  filter  aliasing. 

NCH  indicates  number  of  channels  analyzed  (i.e.,  number 
of  channels  with  8  Hz  filters. 


Fit,.  VI. 3.1  Signal  traces  (100  seconds)  for  Event  4  in  Table  VI. 3.1, 

recorded  at  channel  01A00  with  8  Ha  filter  and  20  Hz  sampling 
rate.  The  unfiltered  data  are  shown  in  the  bottom  trace,  while 
the  first  2  traces  show  the  data  filtered  1.2-3. 2  Hz  and  2.0- 
4.0  Hz,  respectively.  Trace  3  and  4  are  resampled  to  10  Hz 
and  then  filtered  in  the  same  pass  bands.  The  numbers  to 
the  left  of  each  trace  are  maximum  amplitudes,  and  the  traces 
are  scaled  individually. 
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Fig.  VI. 3. 2  SNR  (or  STA/LTA)  curves  for  the  traces  plotted  in  Fig.  VI. 3.1 
(Event  A),  i.e.,  for  2  sampling  rates  (20  and  10  Hz)  and  2 
filters  (1.2-3. 2  and  2. 0-4.0  Hz).  The  numbers  to  the  left 
of  each  trace  show  maximum  SNR,  and  the  traces  are  plotted 
to  the  same  scale. 


filtered  1.2-3. 2  Hz  and  2. 0-4.0  Hz 
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VI. 4  Seismic  noise  at  high  frequencies 

Our  program  for  noise  measurements  in  and  around  the  NORSAR  siting  area, 
reported  upon  also  in  previous  Semiannual  Reports,  has  continued  in 
the  following  way: 

-  All  NORSAR  subarrays  have  been  visited  with  Kinemetrics  PDR-2 
equipment,  'tapping'  the  data  between  the  amplifier  and  the  anti¬ 
aliasing  filter.  The  sampling  rate  mostly  used  has  been  62.5  Hz 
with  filter  at  25  Hz,  but  some  data  have  also  been  recorded  with 
a  125  Hz  rate  and  filter  at  50  Hz.  This  system  has  an  excellent 
dynamic  resolution  up  to  at  least  40  Hz. 

-  The  present  NORESS  area  (subarray  06C)  has  been  studied  in  greater 
detail  (PDR-2  recordings),  using  both  standard  06C  sensors  and  the 
temporary  NORESS  sites.  Time-of-day  variations  have  also  been  studied 
here.  Excellent  dynamic  resolution. 

-  Data  transmitted  in  analog  form  to  NORSAR  from  06C  and  02B  have  also 
been  analyzed,  the  sampling  rate  here  Is  40  Hz,  but  the  dynamic  resolu¬ 
tion  is  poor,  and  there  is  an  Increasing  amount  of  system  noise  from 
the  telephone  lines  for  frequencies  above  2-3  Hz. 

-  Noise  data  from  a  60  m  hole  have  been  recorded  simultaneously  with 
surface  data  (using  the  PDR-2),  and  comparisons  are  made  under  various 
conditions.  Excellent  dynamic  resolution. 

-  Six  other  sites  in  southeastern  Norway  have  been  studied  by  'tapping' 
(using  the  PDR-2)  the  analog  (telephone)  lines  from  the  Southern  Norway 
Seismic  Network.  (SNSN).  System  noise  is  a  problem  also  here. 

Representative  noise  spectra  from  6  NORSAR  subarrays  (except  06C)  are 
shown  in  Figures  VI. 4. 1-2,  with  recordings  using  the  Kinemetrics  PDR-2 
on  the  output  from  the  NORSAR  RA-5  amplifiers  in  the  Central  Terminal 
Vault  (CTV)  of  each  subarray.  This  corresponds  to  System  No.  12  in  Table 
III. 4.1  of  this  report,  where  it  is  seen  that  the  PDR-2  gives  a  gain 
of  42  dB  (for  the  weakest  signals)  in  addition  to  the  72  dB  from  the 


NORSAR  amplifier.  The  recordings  have  been  done  with  a  sampling  rate 
of  62.5  Hz  (16  ms  sampling  interval)  and  with  an  anti-aliasing  filter 
(12  dB/oct)  at  25  Hz.  The  spectra  in  Figs.  VI. 4. 1-2  cover  various  parts 
of  the  winter  from  day  337/1982  to  day  56/1983  and  should  therefore  be 
expected  to  cover  generally  high  noise  levels  for  frequencies  below  2  Hz. 
The  1  Hz  levels  for  the  cases  shown  is  between  5  and  10  dB  above  1  nm^/Hz, 
while  we  previously  have  shown  that  the  typical  summer  level  is  more  in 
the  range  0  to  7  dB.  The  spectra  in  Figs.  VI. 4. 1-2  show  typical  noise 
levels  for  each  subarray,  and  it  is  easily  seen  that  they  are  strikingly 
similar  (above  2  Hz),  with  a  10  Hz  level  at  -50  dB  below  1  nm^/Hz.  The 
only  exceptional  feature  in  these  spectra  is  a  strong  12.5  Hz  noise 
peak  from  02C00,  that  one  is  caused  by  inductive  Interference  from  a 
50  Hz  power  line  close  to  the  analog  tranmlsslon  line  from  seismometer 
to  CTV. 

From  Figs.  VI. 4. 1-2  we  can  also  see  for  most  of  the  channels  some  spectral 
flattening  above  15-20  Hz.  This  is  caused  by  various  cultural  activities 
at  close  distances,  and  is  quite  often  much  more  prominent  than  shown 
here,  sometimes  also  starting  at  frequencies  as  low  as  5  Hz.  There  is  a 
clear  correlation  between  these  observations  and  the  population  density 
in  the  vicinity  of  each  of  the  subarrays. 

Data  recorded  with  a  sampling  rate  of  125  Hz  and  an  anti-aliasing  filter 
at  50  Hz  have  also  been  analyzed,  as  shown  in  Fig.  VI. 4. 3  for  subarrays 
02B  and  06C  (NORESS).  It  is  seen  there  that  the  ambient  noise  level 
continues  to  drop  at  about  the  same  rate  all  the  way  up  to  50  Hz  (pro¬ 
vided  that  cultural  noise  sources  can  be  avoided),  a  rate  which  is  close 
to  15  dB/octave  or  50  dB/decade.  It  is  seen  that  the  lowest  noise 
level  resolved  here  is  about  -80  dB  relative  to  nm^/Hz,  or  260  dB  below 
1  nm^/Hz,  which  is  5  dB  below  Herrin's  (1982)  lowest  value  for  Lajltas. 

The  spectral  differences  between  the  two  sensors  in  Fig.  VI. 4. 3  are 
typical,  in  the  sense  that  subarray  02B  is  somewhat  less  affected  by 
cultural  noise  at  high  frequencies  than  06C,  which  is  the  present  NORESS 
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Because  of  the  possibility  that  the  06C  area  could  be  chosen  as  the  site 
also  for  the  permanent  NORESS  installation,  fairly  extensive  noise  sur¬ 
veys  have  been  conducted  there.  As  a  part  of  that,  we  have  covered  the 
area  with  regular  PDR-2  recordings  throughout  the  day  and  night,  and 
Fig.  VI. 4. 4  shows  here  an  example  for  1400,  1800  and  2200  GMT.  In  addi¬ 
tion  to  the  general  decrease  in  the  noise  level  above  5  Hz  by  night, 
there  is  also  a  prominent  and  fairly  broad  peak  around  5. 5-6.0  Hz  which 
also  disappears  by  night.  This  peak  can  also  be  seen  in  several  of  the 
spectra  presented  above,  and  we  have  come  to  the  conclusion  that  it  prob¬ 
ably  is  caused  by  the  integrated  effect  of  all  the  traffic  in  the  general 
area.  The  survey  has  moreover  shown  that  the  06C  (NORESS)  area  is  at  times 
significantly  affected  by  local  cultural  noise  sources  above  5-10  Hz, 
and  much  of  this  comes  from  lumbering  and  other  economic  activities 
within  the  immediate  vicinity  (<  5  km)  of  the  sites.  In  this  sense,  both 
subarray  02B  and  03C  (which  both  could  be  acceptable  as  NORESS  sites 
from  a  signal  focusing  point  of  view)  are  somewhat  better  than  06C, 
but  the  difference  is  not  very  clear  and  all  seven  subarray  sites  are 
at  times  significantly  affected  by  high-frequency  cultural  noise 
from  sources  at  short  distances. 

All  the  seismic  noise  spectra  presented  above  have  been  taken  from  sur¬ 
face  installations  (i.e.,  2-5  m  borehole  depths).  In  order  to  test  the 
possible  gain  from  deeper  borehole  installations  we  have  conducted  some 
simultaneous  PDR-2  recordings  from  surface  and  60  m  deep  seismometers 
at  site  01A01,  with  results  as  shown  in  Fig.  VI. 4. 5.  It  is  seen  there 
that  the  difference  starts  at  around  8  Hz  and  Increases  to  more  than 
10  dB  for  frequencies  above  20  Hz.  However,  in  comparing  with  the  02B 
and  06C  spectra  in  Fig.  VI. 4. 3  we  see  that  those  spectra  are  closer 
to  the  01A  60  m  spectra  in  Fig.  VI. 4. 5,  and  we  should  therefore  not 
expect  a  gain  of  this  size  for  sites  that  are  less  affected  by  local 
cultural  noise  than  01A.  For  technical  reasons,  we  have  so  far  not  been 
able  to  investigate  the  noise  reduction  potential  for  other  depths  than 
60  m. 
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Seismic  noise  spectra  from  3-component  recordings  done  at  subarray  06C 
are  shown  in  Fig*  VI. 4. 6,  using  SS-1  seismometers  and  the  usual  PDR-2 
recorder.  The  results  show  that  there  is  no  systematic  difference  between 
vertical  and  horizontal  components  for  frequencies  above  3-4  Hz,  while 
for  lower  frequencies  the  noise  level  on  the  vertical  component  is 
slightly  higher. 

We  have,  in  addition  to  the  detailed  noise  studies  from  various  NORSAR 
installation,  also  analyzed  data  from  the  Southern  Norway  Seismic  Network 
(SNSN),  in  this  case  also  using  the  PDR-2  recorder  in  combination  with 
significantly  preamplified  signals.  The  results  clearly  confirm  our 
conclusions  from  above  with  respect  to  the  stability  of  the  noise  spectra 
for  higher  frequencies.  However,  the  data  are  somewhat  affected  by  trans¬ 
mission  noise  on  the  telephone  lines  used,  a  problem  which  also  affects 
our  present  analog  data  channels  from  subarrays  02B  and  06C,  for  fre¬ 
quencies  above  2-3  Hz. 

In  order  to  facilitate  the  comparison  between  a  typical  NORSAR  (or 
southeastern  Norway)  noise  spectrum  and  previously  published  quiet 
sites,  we  have  in  Fig.  VI. 4. 7  plotted  our  resultB  as  presented  above 
on  top  of  the  Queen  Creek  and  Lajitas  noise  spectra  (Herrin,  1982). 

We  see  from  that  figure  that  there  is  not  much  difference  between 
these  spectra  for  higher  frequencies,  and  this  makes  it  natural  to 
ask  the  question  if  there  actually  is  a  fairly  stable  and  uniform 
ambient  noise  level  globally  for  these  frequencies.  In  fact,  since 
southeastern  Norway  and  lajitas  (Texas)  are  so  similar  (as  shown  in 
Fig.  VI. 4. 7),  there  are  no  reasons  to  believe  that  other  areas  of  the 
world  should  be  significantly  different  (provided,  of  course,  that 
local  cultural  noise  sources  can  be  avoided). 

In  conclusion,  we  have  found  that: 

1)  The  ambient  seismic  noise  level  above  2-3  Hz  in  southeastern 

Norway  is  very  stable  both  in  time  and  space.  The  level  at  10  Hz 
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Is  always  very  close  to  -50  dB  relative  to  1  nm^/Hz,  with  a  slope 
in  the  5  to  AO  Hz  range  of  15  dB/octave  (50  dB/decade).  Levels 
down  to  10“®  nm^/Hz  (-80  dB)  have  been  resolved. 

2)  There  is,  for  all  of  the  sites  studied,  a  problem  with  cultural 
noise  for  frequencies  above  5  Hz,  and  this  problem  Increases  with 
increasing  frequency.  There  is  a  rush-hour  peak  around  5. 5-6.0  Hz, 
while  the  cultural  noise  above  10  Hz  sometimes  consists  of  a  broad¬ 
band  contribution  (many  sources  at  various  distances)  but  more 
often  of  narrow  spectral  lines  (few  and  close-in  sources).  This 
problem,  however,  is  manageable  if  careful  site  surveys  are  executed. 

3)  Another  prominent  cultural  noise  source  consists  of  inductive 
interference  from  nearby  50  Hz  power  lines,  and  thiB  often  gives 
significant  contributions,  especially  at  12.5  Hz.  We  know  that 
most  of  this  interference  comes  from  the  cables,  but  we  cannot 
exclude  the  possibility  of  some  interference  also  in  the  seis¬ 
mometer  Itself  if  the  distance  to  the  power  line  is  small. 

4)  In  comparing  the  seven  subarrays  at  NORSAR,  we  find  that  this  50  Hz 
(12.5  Hz)  noise  problem  is  particularly  serious  for  subarrays  01B 
and  02C.  Subarray  01A  is  problematic  with  respect  to  other  cultural 
noise  sources  and  to  some  extent  also  06C  (the  present  NORESS  site), 
while  subarrays  02B,  03C  and  04C  are  somewhat  better  in  this  respect, 
even  though  these  too  are  somewhat  affected  by  the  same  problem.  This 
relative  rating  is  consistent  with  what  we  should  expect  from  merely 
looking  at  the  population  density  and  the  network  of  roads  in  the 
area. 

5)  The  local  cultural  noise  for  frequencies  above  8  Hz  can  be  reduced 
considerably  by  using  60  m  deep  boreholes,  but  we  do  not  know  the 
noise  reduction  potential  for  other  depths.  Wind  noise  is  a  small 
problem  in  this  respect,  partly  because  there  usually  is  very  little 
wind  in  the  area,  and  partly  because  most  of  the  wind  noise  can  be 
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avoided  simply  by  cutting  trees  near  the  site  and  by  Installing  the 
seismometer  at  a  few  meters  depth  in  competent  (crystalline)  rock.8. 

6)  For  higher  frequencies  there  is  no  systematic  difference  between 
vertical  and  horizontal  components. 

7)  In  comparing  with  so-called  'extremely  quiet'  sites  in  other  con¬ 
tinents  (notably  Lajitas  in  Texas)  we  find  that  the  noise  levels 
there  for  frequencies  above  about  10  Hz  are  quite  close  to  what 
we  find  for  southeastern  Norway.  The  possibility  therefore  exists 
that  the  ambient  noise  level  that  we  have  found  for  southeastern 
Norway,  with  its  stability  in  time  and  space,  in  fact  could  be 

a  globally  representative  noise  level. 


H.  B ungum 


Reference 


Herrin,  E.T.  (1982):  The  resolution  of  seismic  instruments  used  in 
treaty  verification  research.  Bull.  Seism.  Soc.  Am.  72,  S61-S67. 


Fig.  VI. 4.1  Noise  power  spectra  (PDR-2  recording  at  62.5  Hz)  for  NORSAR  site  01A05 

(left,  day  55/1983),  01B04  (center,  day  337/1982),  and  02B00  (right,  day  56/1983). 
All  are  day-tine  spectra. 
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FI*;.  VI. 4. 3  Noise  power  spectra  (PDR-2  recording  at  125  Hz)  for  NORSAK 

site  02B00  (left,  day  56/1983)  and  06C02  (right,  day  53/1983) 
Both  are  day-time  spectra. 


61 


Fig.  VI.4.4  Noise  power  spectra  (PDR-2  recording  at  62.5  Hz)  for  NORESS  Site  8  (close  to 
06C02 )  from  day  83/1983  at  1400,  1800  and  2200  GMT. 
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Fig.  VI. 4. 5  Noise  power  spectra  (PDR-?  recording  at  125  Hz)  for  NORSAR 
Site  01A0I  and  from  a  60  m  borehole  at  the  same  site,  at 
two  different  times  (day  55/1983  and  day  87/1983). 


POWER  DENSITY  (NM«*2/HZ) 


Fig.  VI. 4. 6  Noise  power  spectra  (PDR-2  recording)  for  NORSAR  Site 

06C02,  using  3-component  SS-1  seismometers,  day  53/1983. 
Both  are  day-time  spectra,  left :  1314  GMT,  125  Hz  recordin 
ri^ht  :  1349  GMT,  62.5  Hz  recording. 
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VI. 5  Power  spectral  bias  sources  and  quantization  levels 
In  seismic  data  acquisition  systems  there  are  three  main  source-  of  non- 
seismlc  noise  (system  noise)  causing  various  kinds  of  limitations  in 
resolution  and  dynamic  range: 

1.  The  seismometer/amplifier  system 

2.  The  transmission  system  (modulators/demodulators,  telemetry 
equipment,  telephone  lines,  etc.) 

3.  The  recording  system,  for  digital  systems  mainly  quantization 
errors. 

A  typical  problem  associated  with  the  Type  1  noise  source  is  given  a 
seismometer/amplifier  system  with  a  certain  dynamic  range,  to  find  a  gain 
level  that  gives  a  reasonable  tradeoff  between  the  conflicting  needs 
for  resolving  the  seismic  background  noise  and  at  the  same  time  avoid¬ 
ing  clipping  of  strong  signals.  This  is  a  particular  problem  with  regional 
and  local  earthquake  data  as  the  necessary  dynamic  range  for  these  sig¬ 
nals  is  significantly  larger  than  for  teleseismic  events. 

We  have  at  times  experienced  significant  problems  associated  with  the 
Type  2  noise  source  when  using  leased  telephone  lines  for  data  trans¬ 
mission  on  analog  form.  This  applies  in  particular  to  systems  9,  13 
and  14  described  in  Section  III. 2,  where  we  have  found  (for  the  par¬ 
ticular  amplifier  gain  used  there)  that  the  system  noise  is  approxi¬ 
mately  white  and  that  it  reaches  the  level  of  the  ambient  background 
noise  at  a  frequency  of  4-5  Hz,  causing  an  increasing  amount  of 
spectral  bias  for  higher  frequencies. 

For  many  seismic  systems  the  main  limiting  factor  with  respect  to  dynamic 
resolution  is  the  limitations  within  the  recording  system  (Type  3  errors) 
more  than  those  of  the  seismometer/amplifier  system.  For  digital  systems, 
the  main  limitation  is  the  number  of  bits  available  for  representation 
of  signal  level.  This  problem  is  sometimes  'solved'  by  gain-ranging 


(increasing  the  quantization  step  with  increasing  signal  level),  which 
works  well  as  long  as  the  signal  spectrum  is  reasonably  white  and/or 
band-limited  in  frequency.  However,  the  real  situation  with  seismic  sig¬ 
nals  (including  earth  noise)  is  that  weak  high-frequent  signals  are 
superimposed  on  strong  low-frequent  ones,  and  quantization  errors 
are  therefore  dependent  upon  the  shape  of  the  noise  and  signal  spectra 
(for  a  discussion  of  quantization  errors  under  simpler  conditions, 
see  Oppenheim  and  Schafer,  1975). 

We  have  addressed  the  problem  of  quantization  errors  in  a  strictly 
empirical  way,  starting  with  the  unbiased  earth  noise  spectra  presented 
in  Section  VI. 4.  Three  different  power  spectra  were  selected,  where  the 
main  difference  was  that  they  were  based  on  data  sampled  at  125,  62.5 
and  40  Hz,  with  filters  at  50,  25  and  12.5  Hz,  respectively  (correspond¬ 
ing  to  systems  16,  12  and  17  in  Section  III. 4).  Each  of  the  (integer) 
time  series  were  then  scaled  down  successively  by  factors  of  two,  and 
the  power  spectra  calculated.  The  results  are  given  in  Fig.  VI. 5.1, 
where  four  power  spectra  are  calculated; 

1)  Scaling  factor  2°,  with  0.000334  NM/QU  at  1  Hz 

2)  Scaling  factor  28,  with  0.0855  NM/QU  at  1  Hz 

3)  Scaling  factor  210,  with  0.342  NM/QU  at  1  Hz 

4)  Scaling  factor  2^,  with  1.37  NM/QU  at  1  Hz 

It  is  seen  from  the  figure  that  the  power  spectra  for  1)  and  2)  overlap 
almost  completely  (all  possible  effects  of  the  PDR-2  gain-ranging  have 
been  removed  by  a  scaling  of  2^),  at  quantization  level  3)  there  is  a  clear 
bias  extending  in  the  worst  case  almost  down  to  10  Hz,  while  for  level  4) 
the  bias  extends  to  about  5  Hz.  For  these  and  similar  curves  we  have 
derived  the  relationship  shown  in  Fig.  VI. 5. 2,  which  shows  the  minimum 
quantization  level  that  is  required  for  resolving  earth  noise  at  a  par¬ 
ticular  frequency,  provided  that  the  anti-aliasing  filter  is  above  the 
frequency  considered.  It  must  be  emphasized  here  that  this  relationship 
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is  dependent  upon  the  earth  noise  power  spectrum  as  well  as  upon  the 
system  response  function.  If  used  as  a  guideline  for  choosing  quantization 
levels,  it  would  be  wise  of  course  to  select  a  level  a  factor  of  two 
better  than  the  minimum  requirement. 

The  results  presented  above  can  be  used  also  for  investigating  the 
limits  in  resolution  for  the  standard  NORSAR  system,  which  (see  also 
Section  III. A)  have  the  following  very  severe  gain-ranging: 


1) 

Sample 

range 

1-127, 

quantization  0.0427 

NM/QU 

2) 

Sample 

range 

128-511, 

quantization  0.171 

NM/QU 

3) 

Sample 

range 

512-2047, 

quantization  0.683 

NM/QU 

4) 

Sample 

range 

2048-8191, 

quantization  2.73 

NM/QU 

We  find  at  NORSAR  that  the  peak  amplitudes  of  the  seismic  background  noise 
are  usually  at  sample  range  2)  but  occasionally  at  range  3),  which  means 
that  the  quantization  is  either  4  or  16  times  poorer  than  indicated  by 
the  usually  quoted  value  of  0.0427  NM/QU.  In  order  to  test  this  on  real 
data  we  have  in  Fig.  VI. 5. 3  computed  three  noise  power  spectra  covering 
the  same  time  interval,  with  1)  unbiased  PDR-2  data,  2)  NORSAR  data 
filtered  at  8  Hz,  and  3)  NORSAR  data  filtered  at  4.75  Hz.  These  filters 
are  very  sharp  (24  dB/octave),  and  we  can  never  recover  unbiased  data 
above  the  cutoff  frequency.  In  fact,  from  Fig.  VI. 5. 3  we  see  that  the 
bias  in  the  worst  case  may  extend  down  to  between  3  and  4  Hz.  However, 
the  data  used  in  Fig.  VI. 5. 3  are  taken  from  different  channels  causing 
some  level  differences  also  for  lower  frequencies,  and  it  is  therefore 
not  possible  to  find  a  more  exact  'bias  frequency'  from  such  tests 
(there  are  also  some  instabilities  in  time). 

Another  implication  of  the  results  presented  here  is  that  an  even  stronger 
bias  should  be  expected  for  earthquake  and  explosion  signals  with  peak 
amplitudes  in  sample  range  3)  and  a  sharp  rolloff  in  the  spectrum  towards 
higher  frequencies. 
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The  best  way  to  solve  these  problems  would  of  course  be  to  replace  the 
AD-converter  with  a  more  advanced  one  (Increasing  significantly  the 
number  of  bits),  but  some  improvement  (for  small  signals)  within  the 
frame  of  the  old  system  could  also  be  obtained  by  whitening  the 
uncorrected  power  spectrum  through  the  introduction  of  an  extra 
analog  high-pass  filter  (with  a  reasonably  gentle  slope). 

H.  Bungum 

Reference 

Oppenheim,  A.V.  and  R.V.  Schafer  (1975):  Digital  Signal  Processing 
(Chapter  9),  Prentice-Hall,  New  Jersey. 
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Quantization  (NM/QU) 


Fig.  VI. 5. 2  Relationship  showing  the  frequency  range  for  which  suf¬ 
ficient  noise  resolution  can  be  obtained  (unbiased  spectra), 
given  a  certain  quantization  level  at  1  Hz.  The  relation¬ 
ship  is  valid  only  for  the  particular  combination  of  noise 
power  spectrum  observed  in  southeastern  Norway  (see 
Section  VI. 4)  and  the  response  functions  of  the  systems 
used  to  delineate  it  (see  Section  III. 4,  systems  12,  16 
and  17). 


POWER  OENSITY  (NM*»2/HZ) 
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Fig.  VI.  5. 3 


Noise  power  spectra  for  day  62/1983  (left),  and  day 
87/1983  (right),  and  with  each  case  covered  by:  1) 
an  unbiased  PDR-2  recording  with  sufficient  quantiza¬ 
tion  (62.5  Hz  sampling  rate,  25  Hz  filter);  2)  NORSAR 
20  Hz  recording  with  8  Hz  filter  (center  trace);  and 
3)  NORSAR  20  Hz  recording  with  5  Hz  filter  (top  trace, 
clearly  biased). 
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VI. 6.  Magnitudes  from  P  coda  and  Lg  using  NORSAR  data 

The  objective  of  this  study  is  to  investigate  the  stability  of  magnitude 
estimates  based  on  P  coda  and  Lg,  compared  to  that  of  conventional  P- 
wave  magnitudes.  Using  NORSAR  recordings  of  presumed  explosions  from 
Semipalatinsk  as  a  data  base,  the  purpose  has  been  in  particular  to 
evaluate  amplitude  variations  across  the  array  as  well  as  compare  NORSAR 
magnitude  estimates  to  those  of  ISC  and  NEIS. 

For  the  purpose  of  this  study,  NORSAR  recordings  of  25  presumed  explosions 
from  Semipalatinsk  (Table  VI. 6.1)  were  used  as  a  data  base.  The  epicentral 
distance  is  about  38  degrees  for  these  events,  and  Lg  is  usually  only 
slightly  stronger  than  the  preceding  P  coda  (Fig.  VI. 6.1).  It  is  noteworthy 
that  the  coda  energy  stays  significantly  above  the  noise  level  for  several 
minutes  following  Lg,  in  the  case  of  high  magnitude  events.  However,  we 
have  found  that  below  mg  «  5.5  the  noise  influence  on  the  Lg  arrival  be¬ 
comes  significant  except  during  very  quiet  noise  conditions. 

We  will  first  consider  the  amplitude  variations  across  NORSAR  for  a  given 
event.  The  significant  variability  of  P  amplitudes  across  NORSAR  is  well 
known,  as  illustrated  in  Fig.  VI. 6. 2.  The  standard  deviation  is  frequency 
dependent,  and  typically  ranges  from  about  0.20  to  0.25  mg  units.  Further¬ 
more,  the  amplitude  pattern  is  strongly  regionally  dependent  (Berteussen 
and  Husebye,  1974;  Ringdal,  1977).  Fig.  VI. 6. 3  shows  that  the  variability 
of  the  P  coda  is  significantly  lower  (between  0.07  and  0.10  mg  units 
standard  deviation  of  peak  amplitudes  in  30  second  windows  using  a  0. 6-3.0 
Hz  filter)  and  in  fact  remains  fairly  constant  out  to  and  beyond  the  Lg 
arrival.  Furthermore,  the  variability  of  the  P  coda  across  NORSAR  is 
similar  to  that  of  the  noise  preceding  P,  if  equal  time  windows  are 
used.  If  RMS  values  are  considered,  the  standard  deviations  are  slightly 
lower,  and  average  about  0.06  mg  units. 

We  may  therefore  conclude  that  considering  P  coda  or  Lg  reduces  the  receiver 
'focusing'  effects  in  the  sense  that  the  standard  deviations  across  NORSAR 
are  reduced  by  a  factor  between  2  and  3  compared  to  the  first  few  cycles 
of  P. 
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The  near  receiver  'focusing'  effects  observed  for  P-waves  across  NORSAR 
might  be  expected  to  have  counterparts  in  near-source  'focusing'.  Although 
the  actual  physical  interpretation  is  uncertain,  it  is  nevertheless 
clear  that  significant  regionally  dependent  bias  values  are  observed  at 
NORSAR  even  within  a  very  limited  source  area.  Fig.  VI. 6. 4  shows  NORSAR  mg 
plotted  against  ISC  mg  (or  NEIS  mg  where  ISC  data  are  not  available). 

The  most  noteworthy  feature  is  the  consistently  low  NORSAR  magnitudes 
for  Degelen  Mountains  events  compared  to  events  from  Shagan  River.  However, 
there  is  considerable  variability  even  within  the  latter  region.  For  the 
entire  event  set,  the  NORSAR  bias  (at  instrument  01A04)  varies  from  0.0 
to  about  0.7  mg  units,  with  a  mean  of  0.42  and  a  standard  deviation  of 
0.21.  Since  the  amplitude  pattern  across  NORSAR  is  fairly  stable  for  Semi- 
palatinsk  events,  these  relative  bias  values  remain  also  when  averaging 
magnitudes  across  the  NORSAR  array. 

Fig.  VI. 6. 5  is  similar  to  Fig.  VI. 6. 4,  except  that  NORSAR  mg  is  replaced  by 
log  RMS  Lg  amplitudes  at  NORSAR.  The  standard  deviation  is  reduced  to 
0.10  logarithmic  units,  i.e.  ,  about  half  of  that  obtained  using  conven¬ 
tional  mg.  Furthermore,  the  Degelen  Mountain  events  now  fall  into  the 
general  trend  of  Shagan  River  events.  As  shown  in  Figure  6,  the  scatter 
between  average  and  single  instrument  Lg  observations  is  small. 

We  can  conclude  that  magnitudes  based  on  Lg  effectively  reduce  bias  due 
to  near-receiver  as  well  as  near-source  effects.  The  question  of  how 
accurately  mg  estimates  can  be  obtained  using  Lg  from  a  single  station 
(or  an  array)  for  a  region  like  Semlpalatinsk  remains  open,  but  it  seems 
reasonable  that  the  standard  deviation  of  0.10  mg  units  represents  aji 
upper  bound.  Indeed,  it  is  quite  common  that  mg  estimates  by  NEIS  and  ISC 
differ  by  0.1  mg  units  or  more,  even  though  many  of  the  same  stations  are 
used.  Furthermore,  about  half  of  the  stations  used  in  these  estimates  are 
located  in  Central  Europe,  thus  implying  that  the  NEIS  and  ISC  estimates 
may  be  severely  biased  because  of  the  lack  of  homogeneous  geographical 
station  distribution.  Thus  it  is  quite  possible  that  single  station  mg 
(Lg)  magnitudes  may  have  an  accuracy  at  least  as  good  as  those  of  NEIS 
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or  ISC.  Additional  Improvements  may  of  course  be  obtained  by  averaging 
Lg  magnitudes  from  a  well-distributed  network,  of  stations. 

F.  Ringdal 
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2/16/73  5  2  57  7  49.835N  78.232E  0 

8-BP  2.0-4. 0  10HZ  3RD 


Fig.  VI. 6. 2  Typical  P-vave  amplitude  distribution  across  NORSAR  for 
Semipalatinak  events.  The  figure  covers  30  seconds  of 
filtered  SP  data  (2-4  Hz)  for  each  center  sensor  of  the  22 
NORSAR  subarrays.  Amplitudes  vary  by  a  factor  of  10  for  the 
initial  P  onset,  whereas  the  amplitude  variation  Is  less 
pronounced  in  the  P  coda. 


NORSAR  P  CODA  DECAY  -  SHAGAN  RIVER 


VI. 6. 3  Typical  P  coda  decay  of  NORSAR  recordings  from  Shagan  River 
events*  The  standard  deviations  of  peak  amplitudes  across 
NORSAR  in  30  second  windows  are  indicated  by  vertical  bars. 
Data  for  the  plot  has  been  obtained  by  combining  Events  2 
and  8  of  Table  VI. 6.1,  and  a  filter  of  0.6-3.0  Hz  has  been 
applied. 
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N-20  BIAS-  0.42  STD-0.206 


MAGNITUDE  (MB)  -  ISC.NEIS 


Fi^.  VI. 6. 4  NORSAR  mg  (instrument  01A04)  versus  ISC  or  NEIS  mg  (ISC  data 
used  when  available).  Note  the  difference  in  m^,  bias  between 
events  from  Shagan  River  and  Degelen  Mountains.  The  stippled 
lines  correspond  to  plus  and  minus  two  standard  deviations. 


N-25  BIAS-  2.94  STD-0.099 


Fig.  VI. 6. 5  Same  as  Fig.  VI. 6. 4,  except  that  NORSAR  m^  has  been  replaced 
by  an  Lg  energy  estimate.  This  estimate  has  been  obtained  by 
averaging  log  RMS  of  individual  NORSAR  sensor  traces  filtered 
in  the  band  0. 6-3.0  Hz.  The  time  window  is  2  minutes,  starting 
40  seconds  before  Lg  arrival.  Note  that  Shagan  River  and 
Degelen  Mountain  events  now  show  good  consistency. 


m 


-  80  - 


> 


o 


— J 

pH 

»4 

m  ON 

m 

u>  o  in  cn 

omc'f'*'rcn<tf’cvp**u> 

m  O' 

m  cm  <)  ui 

<-> 

r- 

CM  O 

O' 

MJ 

CM 

go  r-»  uj 

M' 

ooO'O'ooo'nir— 

00 

O'  cu 

WtTMJ 

u> 

n> 

• 

• 

• 

• 

• 

• 

•  •  • 

• 

•  •  •  • 

•  «  • 

« 

• 

•  » 

•  *  • 

• 

O 

O 

CM 

»M 

CM 

CM 

CM 

<M 

2 

2 

2 

CM 

2 

2 

2 

3 

cn  cu  cm 

CM 

CU  CM  CM  m  CM  CM  m 

— i 

or* 

m  o  r* 

r»i 

tu 

M 

4 

uj  4  r*  cn 

4  r-  co 

Ui 

c- 

O#  Cl 

ui  us  r- 

O 

4* 

4 1 

r- 

in 

Ml 

at 

uj  tn  ui 

Ui 

cn  cn  Ui  u> 

U3  U>  4 

4  u> 

cn  u- 

p*  Ultfim 

r- 

w  w  o 

0 

c 

c 

a 

a 

c 

0 

0 

0 

c 

c 

G 

0 

C 

c 

0 

0 

c 

c 

0 

c 

a 

m 

• 

• 

• 

• 

• 

• 

•  •  • 

• 

•  *  •  * 

#  «  • 

• 

• 

•  • 

•  •  • 

• 

u 

u  a 

o  v»u 

O  O  W 

c 

0 

0 

c 

c 

c 

G 

C 

c 

0 

0 

c 

c 

0 

c 

c 

o 

> 

_J 

< 

u 

u  cu 

in 

oj 

sU 

UP  Ul  O' 

00 

O'  CM  CM  Ul  O  CT  CO 

in  O'  o  w 

u-  u*  a 

cc 

Ui 

pi  cr 

O' 

PI 

CD  00  41 

O'  CO  -«  o 

h  tr  c*i 

r* 

O'  o  tr 

«  o  w 
•  •  • 

o 

cc 

• 

• 

• 

• 

• 

• 

•  •  • 

• 

•  •  •  • 

•  •  • 

• 

• 

•  • 

• 

C3 

< 

*M 

14 

CM  CM 

CM 

*M  CM  CM  CM 

CM 

cm  m  m  pi 

cn  CM  CM 

«M 

CM 

cn  cm 

ni  cn  cm 

Cl 

-1 

a: 

>*• 

>r 

CM  O 

4 

CJ  CO  4  4 

UJ 

O  CD  O'  CM  CM  CM  r— 

O 

cc  *-• 

<. 

ni 

Ul 

<r  o 

r- 

Ul 

n  n  r*- 

4 

m  u>  m  r- 

4  cn  ui 

CM 

cn  ui 

ad 

i/> 

• 

* 

• 

• 

■ 

• 

•  •  • 

• 

•  •  •  • 

•  •  • 

1 

« 

«  • 

1  1  t 

1 

X 

a: 

MJ 

Ul 

43 

vO 

Ul 

Ul 

6 

6 

c 

UJ 

6 

6 

C 

6 

6 

6 

c 

U> 

UJ  UJ 

a 

z 

o 

a> 

n» 

00  O' 

*C3 

m 

9 

B 

7 

0 

6 

2 

8 

1 

o  a  «"> 

r>* 

aj 

on 

• 

• 

• 

% 

• 

• 

•  »  « 

• 

*  *  •  » 

•  •  • 

• 

| 

1  1 

1  1  1 

1 

X 

•— i 

u  > 

Ul 

in 

Ul 

Ul 

Ul 

in  in  ui 

41 

5 

6 
c 

6 

U)  UJ  ui 

ui 

HJ 

pi 

O' 

O'  in  p> 

9 

8 

7 

0 

6 

3 

8 

1 

o  o  U1 

•  •  • 

r- 

O' 

CMTOOO 

CM 

LU 

M 

• 

• 

• 

• 

• 

• 

•  •  • 

• 

•  •  •  » 

• 

• 

•  * 

•  •  • 

• 

a.  u/  minuunm>nininin<oin'0<n'<M)'OinuMninin«oin'0 
a: 


UJ  vil  Ui  Ui  VU  Ui  UJ  UJ  VU  U4  UJ  UJ  Uy  Uj  UJ  UJ  Ui  Ul  Uj  UJ  Ui  UJ  UJ  tU  UJ 

vDf-HM4,in®<o©mmin‘oo*o^oo'tf\inH-<hN'0 

^  ^('i^^f^MtlO'O^flOH'Om^^O’TCOOOr-MAN 

Q  03  O'  oO^(D(J'HO»(M7'00'0®HCOOO^(MDeD(C 
•••#••••••••••##••••••••• 

WC0OB®®®C0WO®a)0'a30'ffl<DW0'0»a3W®C0W 

422Z4£ZZX«ZZ2^2Z22ZZ22ZZ^ 

44CMfn400'0'4®u3®cn-Mr*«>*o'0®eupuc<Up*r,-cn 

n  a*OOfl5M^«c0h®O5'O^^®N^0'0»»*a9«fl00' 


O' 

o 

o 

O' 

O' 

O' 

O' 

O' 

O' 

O' 

o 

O' 

o 

O' 

O' 

O' 

O' 

O' 

O' 

O' 

O' 

O' 

O' 

O' 

O' 

*• 

in 

in 

'T 

>1* 

4* 

>1* 

* 

4- 

in 

* 

m 

** 

-t 

-4- 

4* 

■i- 

•4* 

4- 

4- 

4* 

4- 

UJ 

Is- 

r* 

r- 

c- 

r'- 

f*- 

r- 

r- 

r* 

h- 

r* 

r- 

'O 

r- 

r- 

r*. 

>r 

>o 

PU 

03 

CM 

® 

4 

X 

tn 

in 

in 

m 

m 

in 

in 

in 

in 

in 

in 

in 

m 

m 

m 

in 

n 

m 

p4 

o 

H 

H 

O 

o 

H 

V- 

u> 

43 

M 

IM 

u> 

43 

43 

43 

43 

4> 

in 

43 

4) 

43 

43 

>0 

43 

ni 

4 

Is- 

n- 

r- 

h- 

in 

cn 

m 

o 

cn 

in 

o 

m 

4* 

4* 

cn 

a 

in 

4 

m 

-H 

cn 

in 

«n 

cn 

4 

pU 

p4 

in 

cn 

4 

Q£ 

cu 

CC 

cc 

<n 

cn 

cn 

cu 

CU 

CU 

cu 

in 

cu 

cn 

cn 

cn 

4 

<n 

CM 

m 

cn 

CM 

cn 

rn 

m 

O 

o 

o 

a 

o 

O 

O 

o 

O 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

O 

O 

O 

o 

r- 

C" 

H 

CO 

00 

00 

00 

00 

00 

CO 

O' 

O' 

O' 

O' 

O' 

o 

o 

o 

o 

p4 

pU 

P«i 

p4 

c- 

r* 

n- 

c- 

c- 

r- 

r* 

f- 

r~ 

c- 

r- 

r- 

f4 

h- 

00 

® 

CD 

00 

00 

CD 

® 

® 

® 

UJ 

*4 

> 

'■s 

*4 

*4 

>» 

4. 

"4 

4. 

*4 

4. 

*4 

4* 

4, 

4* 

X 

V 

4 

4 

4 

4 

4 

4 

c- 

o» 

o* 

in 

in 

UJ 

M 

r-C 

m 

€0 

m 

4 

<n 

Is- 

4 

00 

cu 

cu 

O' 

cu 

4 

cu 

cn 

00 

O' 

h- 

< 

cu 

CM 

n 

cu 

M 

CM 

^4 

o 

cu 

«<u 

O 

ru 

o 

O 

cu 

o 

CM 

CM 

PU 

*u 

CM 

t-4 

CM 

CM 

Q 

V 

'M 

*4 

4, 

4. 

4, 

4* 

V 

4. 

"4 

4. 

*4 

^4 

4. 

4* 

4, 

4* 

•4 

4 

4 

4 

4 

4 

4 

4 

UT 

UJ 

O' 

4* 

cn 

sr 

43 

c- 

Is- 

vn 

4J 

r- 

Oi 

o 

cu 

in 

Ui 

o 

«M 

4 

m 

o 

pU 

CU 

”> 

o 

D 

3 

n 

o 

O 

o 

o 

a 

*-4 

9 

9 

o 

r-c 

pU 

o 

9 

p^ 

C3 

o 

H 

H 

pU 

Z 

UJ 

r-4 

cu 

m 

'T 

in 

UJ 

r- 

a> 

O' 

o 

eg 

cn 

4 

m 

u> 

r* 

OD 

O' 

9 

pU 

CM 

m 

4 

in 

M 

«■* 

pU 

•4 

-4 

cu 

CM 

IM 

CM 

CM 

«M 

UJ 


U 

a> 

OJ 

CO 

-a 

-c 

P-1 

4-1 

■u 

4-1 

•H 

X> 

CO 

c: 

C 

.c 

O0 

<0 

4-1 

cO 

6 

a i 

< 

4-> 

p 

C/3 

o 

a> 

o£ 

z 

0 

O 

ss 

• 

O 

rH 

(0 

o 

c* 

0) 

u 

'T* 

o 

<n 

00 

w 

o 

‘ri 

o 

u 

CO 

4-1 

*0 

03 

p 

03 

E 

(V 

0C 

3 

4-1 

cO 

U 

CO 

V4 

4-1 

TJ 

(33 

(A 

> 

C 

p 

(0 

•H 

• 

0 

cp 

Q 

c 

Jl 

O 

U  00  X? 

C  £  4) 

03  52  Vf  • 

>  3  co 

ai  oo  «  w 

O  (0  c 

0)  f-<  4) 

a>  0  > 

>  *  a» 

r4  ^  Op 

ZQ<S  J  r~» 

O  rH 

Q>  <C  1/3  CO 
-*  -t  i: 
o  o  S 

*  o 

U  *J  CO  4-i 

c  o 

<U  Q)  r-4  a» 

£2  S  -< 
H  3  TJ  J3 
b  C  (0 

•  4J  Ctf  r~H 
>>  «  **"< 
■o  C  C  (0 
3  -H  O  > 

*J  *H  ttf 
CO  w 

«  to  u 

CA  <  i-<  O 
H  W  >  d 
£2  Dd  a» 

*J  O  T)  II 
SB  W 

C  X)  (TJ 

rltj  U 

C  CO  to 
■Q  (0  tj  W 
03  C  3 
to  U  10  H 
3  W  *J  (0 
H  to  > 
(A 

*J  •  500 

C  C/3  C  B 
0)  M  *ri 

>  W  TJ  0$ 

<u  z  c  < 
O  c/1 
44  jn  a.  od 
O  £3  »  O 

o>  as 

*J  ^  u 

to  Ala  33 
■H  0  O  C 

pJ  w  Xj  to 


> 

<0 

H 

Xi 

t2 


VI. 7  A  new  regional  array  In  Norway;  Design  work. 

A  prototype  regional  array  will  be  installed  in  Norway  during  the  summer 
of  1984.  According  to  firm  plans  as  of  1  June  1983,  the  array  will  com¬ 
prise  25  short  period  elements,  4  out  of  which  will  be  3-component 
deployments  with  the  remaining  21  elements  consisting  of  vertical 
motion  seismometers  only.  In  addition,  there  will  be  a  broad-band 
three-component  system.  Sampling  rates  will  be  40  Hz  for  the  short- 
period  channels  and  10  Hz  for  the  broad-band  system. 

The  on-line  processing  of  data  from  the  25  vertical  short-period  chan¬ 
nels  will  be  based  on  an  existing  program  package  (Mykkeltvelt  et  al, 
1982;  Mykkeltvelt  and  Bungum,  1982).  The  subject  of  this  contribution 
is  the  geometrical  configuration  of  the  25-channel  vertical  array. 

Design  experiments 

We  have  previously  devised  a  method  for  array  configuration  optimization 
with  respect  to  SNR  gain  by  beamforming  (Mykkeltvelt  et  al,  1983).  We 
have  demonstrated  optimized  geometries  leading  to  theoretical  gains 
well  in  excess  of  the  standard  /N  gain,  by  utilizing  negative  minima 
in  the  observed  noise  correlation  curves.  Such  optimum  geometries, 
however,  tended  to  be  rather  'peaked'  in  their  frequency  response, 
i.e.,  a  very  high  gain  at  one  particular  frequency  was  generally  ac¬ 
companied  by  low  gains  at  other  (relevant)  frequencies.  The  optimized 
geometries  were  characterized  by  one  particular  intersensor  spacing 
being  represented  as  many  times  as  possible  in  the  geometry.  This  dis¬ 
tance  reflects  the  separation  for  which  the  noise  correlation  has  its 
minimum,  for  a  given  frequency  interval.  For  optimization  taking  several 
frequency  bands  into  consideration  (e.g.,  giving  equal  weight  to  each  of 
five  different  frequency  bands  in  the  gain  expression),  again  one 
single  Intermediate  frequency  'dominated'  the  geometry.  For  on-line 
processing  of  regional  events  on  the  new  array,  signal  frequencies  in 
the  range  1.5  to  5.0  Hz  will  be  of  Importance.  In  this  range  the  dis¬ 
tance  corresponding  to  the  noise  correlation  minimum  varies  by  as  much 
as  a  factor  of  3  and  optimum  geometries  for  different  frequencies 


-  82  - 


within  this  range  would  be  vastly  different.  The  configuration  to  be 
finally  deployed  must  have  many  combinations  of  sensor  pairs  at  optimum 
separation  for  any  frequency  within  a  fairly  wide  range. 

From  the  foregoing  discussion  it  follows  that  we  have  not  been  able  to 
make  much  use  of  our  optimization  procedure  during  the  more  recent 
stages  of  planning  for  the  new  array.  Rather,  a  design  idea  set  forth 
by  Followill  and  Harris  of  1LNL  (Followill  and  Harris,  1983)  has  been 
pursued.  They  propose  a  geometry  based  on  concentric  rings  spaced  at 
log-periodic  intervals  in  radius  R,  according  to  the  formula: 

R  =  ^min  *  *  n  '  0.1.2,...  (1) 

Their  design  includes  the  deployment  of  an  odd  number  of  elements  sym¬ 
metrically  distributed  in  azimuth,  and  it  has  the  following  attractive 
features : 

With  an  odd  number  of  elements  in  each  ring,  the  corresponding 
coarray  (defined  as  the  set  of  all  intersensor  separations,  in 
vector  space)  pattern  has  no  overlap  among  its  points,  i.e.,  it 
samples  the  wavefield  in  the  best  possible  way,  in  this  respect. 
Designs  based  on  (1)  comprise  comprehensive  subsets  of  sensors 
with  very  different  typical  intersensor  separations,  implying 
that  both  high-frequency  and  low-frequency  phases  could  be  well 
enhanced  by  appropriate  subsets  of  the  array. 

-  The  beam  patterns  for  the  above  designs  are  favorable,  with  a 
narrow  main  lobe  yielding  good  resolution  in  phase  velocity  and 
azimuth  and  absence  of  cumbersome  side  lobes. 

More  specifically,  the  configuration  in  Fig.  VI. 7.1,  OR13579,  was 
proposed  for  the  new  array.  It  is  the  realization  of  (1),  with 
Rmin  “  200  m,  a  »  2.25,  n  ■  0,1, 2, 3  and  with  3,5,7  and  9  elements 
in  each  ring,  plus  one  in  the  center.  This  gives  an  array  of  aperture 
about  4.45  km.  For  high-frequency  phases  (3  Hz  and  above)  the  outer 
ring  does  not  contribute  to  the  gain  and  should  be  omitted  during 


processing.  This  leaves  an  array  of  aperture  2  km.  Similarly,  the 
two  inner  rings  should  be  masked  while  processing  low-frequency  phases. 

SNR  gains  by  beamforming  for  the  OR13579  design  and  relevant  sub¬ 
geometries  of  it  were  checked  using  correlation  curves  for  signals 
and  noise  based  on  recordings  on  the  NORESS  array.  Signal  and  noise 
correlations  were  measured  in  the  two  bands  1-3  and  3-5  Hz  and  theo¬ 
retical  gains  were  computed  (Mykkeltveit  et  al,  1983)  for  a  range  of 
values  of  R„iin  in  (1).  For  each  of  the  two  frequency  bands  and  cor¬ 
responding  full  set/subset  of  sensors,  the  proposed  value  of 
(200  m)  gives  gains  close  to  optimum  for  this  design. 

The  beam  pattern  of  OR13579  is  shown  in  Fig.  VI. 7. 2a.  The  narrow 
main  lobe  is  partly  due  to  the  large  aperture  of  4.45  km  of  this 
design.  Standard  beam  pattern  computations,  however,  assume  identical 
signals.  In  order  to  evaluate  the  role  of  the  aperture  in  a  wavenumber 
resolution  context,  we  compute  response  patterns  incorporating  realistic 
signal  correlations,  and  one  example  of  a  resulting  beam  pattern  is 
shown  in  Fig.  VI. 7. 2b.  As  can  be  seen,  the  resolution  capability 
of  the  array  is  degraded  through  the  widening  of  the  main  lobe. 

Because  of  the  above  concerns  about  actual  resolution  capability  for 
a  4.45  km  aperture  array  and  several  logistic  constraints,  the  decision 
was  made  to  plan  for  an  aperture  of  about  3  km  for  the  new  array,  but 
at  the  same  time  retain  the  basic  ideas  underlying  the  OR13579  design. 
The  new  proposed  design,  OR13579  1984,  follows  from  (1)  with  R,„in  * 

150  m  and  a  m  2.15,  yielding  an  outer  ring  diameter  of  2.985  km. 

OR13579  1984  was  checked  with  respect  to  theoretical  gains  in  the 
same  way  as  described  above  for  OR13579.  For  this  task,  we  used  noise 
correlation  curves  for  the  frequency  bands  1.5-2. 5  and  3. 5-4. 5  Hz  in 
conjunction  with  the  previously  derived  signal  correlation  curves  for 
the  bands  1-3  and  3-5  Hz,  thus  making  the  noise  correlation  curves 
more  representative  for  the  center  frequency  in  the  signal  bands. 
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(This  amounts  to  partly  compensating  for  the  strong  decay  with  fre¬ 
quency  In  the  noise  spectra. ) 

Gains  for  the  OR13579  1984  geometry  as  a  function  of  are  given 

in  Fig.  VI. 7. 3  for  different  weighting  of  frequency  bands  and  dif¬ 
ferent  sensor  masking  schemes.  It  is  found  that  for  the  proposed  value 
of  150  m  for  Rmtn  it  does  not  pay  to  delete  the  outer  ring  of  9  elements 
for  the  higher  frequencies.  For  the  lower  frequency  range,  however,  gain 
improvement  is  achieved  by  masking  the  inner  two  rings.  A  general  im¬ 
pression  from  this  figure  is  that  one  should  be  able  to  improve  the 
gain  by  a  slight  increase  of  and  a  corresponding  increase  in  the 

aperture  to  close  to  4  km,  but  the  net  gain  from  this  change  would 
amount  to  less  than  1  dB,  which  is  of  minor  importance  compared  to  the 
above  considerations  on  array  resolution  and  logistic  implications. 

Other  geometricl  patterns,  both  previous  array  realizations  (CPO,  UBO, 
LASA)  and  new  concepts,  have  been  Investigated  along  the  same  lines 
as  above.  None  of  these,  however,  produced  both  a  beam  pattern,  a  co¬ 
array  pattern  and  theoretical  gains  equal  to  or  better  than  the  Followill 
and  Harris  odd-ring  designs. 

1983  temporary  field  Installations 

There  will  be  two  experiments  during  the  summer/fall  of  1983  to  evaluate 
our  current  ideas  of  array  design.  These  experiments  are: 

1)  Temporary  deployment  of  five  three-component  sets.  This  experiment 
is  designed  to  give  us  the  background  data  for  the  decision  of 
where  to  deploy  the  three-component  stations  in  the  1984  array. 

2)  Deployment  of  OR13579  1984  with  the  exception  of  4  instruments 
in  the  outer  ring.  This  gives  us  the  opportunity  of  experimental 
verification  of  the  potentials  of  the  OR13579  1984  design.  There 
is  an  option  for  one  reconfiguration  of  the  21  channels  in  August 
1983. 


S.  Mykkeltveit 
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Flj,.  VI. 7. 2  (a)  Standard  beam  pattern  (impulse  and  response)  for  OR13579 

This  corresponds  to  signal  correlation  equal  to  1  for 
all  intersensor  separations. 


Fig.  VI. 7. 2  (b)  Array  response  for  OR13579  for  realistic  signal  correla¬ 

tions.  Signal  correlation  is  assumed  to  be  linear  and 
equal  to  0  at  5  km  intersensor  separation.  This  value  is 
typical  of  Lg  and  high-frequency  (above  ~  4  Hs)  Pn. 


Gain  (dB)  Gain  (dB)  ^  __  Gain  (dB) 
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Fig.  VI. 7. 3  Theoretical  gains  for  the  OR13579  1984  design  and  relevant 

subgeometries,  as  a  function  of  R^in*  'frequency  weights  1  1* 
means  that  the  low  and  high  frequency  bands  (1-3  and  3-5  Hz, 
respectively)  are  given  equal  weight  in  the  gain  estimation,  a.s.o 


